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ABSTRACT 


The  feasibility  of  utilizing  detonations  for  air-breathing  propulsion  is  the  subject 
of  a  significant  research  effort  headed  by  the  Office  of  Naval  Research.  Pulse  Detonation 
Engines  (PDE)  have  a  theoretically  greater  efficiency  than  current  combustion  cycles. 
However,  pulse  detonation  technology  must  mature  beginning  with  research  in  the 
fundamental  process  of  developing  a  detonation  wave.  This  thesis  explores  various 
ignition  conditions  which  minimize  the  deflagration-to-detonation  transition  distance 
(XDDt)  of  a  single  detonation  wave  in  a  gaseous  mixture. 

Specifically,  the  minimum  XDDT  was  determined  for  different  Ethylene  and 
Oxygen/Nitrogen  gaseous  mixtures  under  varying  ignition  energy  (0.33-8.31  Joules), 
mixture  equivalence  ratios  (0.6-2.0),  and  ignitor  locations.  To  conduct  the  experiments  a 
6ft.  long,  3in.  diameter  tube  combustor,  support  equipment,  and  operating  software  was 
built.  Four  independent  test  scenarios  were  investigated  and  trends  developed  to 
determine  the  minimum  XDDT  while  reducing  oxidizer  blend  ratios. 

Results  show  that  XDDT  significantly  depends  on  mixture  equivalence  ratio  ((p)  and 
was  minimized  at  cp  ~  1.1.  No  dependence  on  ignition  energies  greater  than  0.5  Joules 
was  observed.  A  further  reduction  in  Xddt  was  observed  with  the  ignitor  located  one 
combustor  diameter  from  the  head  wall.  These  results  will  be  useful  in  future  designs  of 
pre-detonators  for  larger  PDEs. 
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I. 


INTRODUCTION 


A.  BACKGROUND 

Almost  all  current  propulsion  systems  utilize  a  constant  pressure  (Cp)  combustion 
process  to  add  energy  to  the  working  fluid.  Most  engine  designs  are  based  on  the 
constant  pressure  burning  of  the  fuel  within  the  combustor.  An  approach,  though  not 
new,  is  emerging  to  investigate  engine  efficiency  limits  within  the  combustion  process 
itself.  This  idea  is  based  on  the  fact  that  fuels  bum  significantly  more  efficiently  in  a 
nearly  constant  volume  compression  based  event  verses  a  constant  pressure  process.  In  a 
practical  flowing  system,  this  is  commonly  referred  to  as  a  detonation.  The  practicality 
of  applying  a  detonation  to  a  propulsion  system  lies  in  the  ability  to  couple  this  increase 
in  thermal  efficiency  to  an  increase  in  propulsive  efficiency.  To  do  this,  the  physics 
behind  establishing  a  detonation  and  its  propagation  along  a  combustor  axis  must  be 
understood  and  harnessed  before  it  can  be  used  effectively  in  a  controlled  manner  for 
propulsion  applications. 

B.  HISTORY  OF  DETONATION  RESEARCH 

1.  Discovery  of  the  Detonation 

The  first  recognized  detonation,  then  called  a  ‘violent  explosion’,  was  discovered 
and  later  patented  by  A.  Nobel  in  1863.  He,  with  his  father,  invented  a  mercury 
fulminate  ignitor  that  initiated  a  detonation  in  a  nitroglycerine  charge.  He  later  perfected 
this  process  in  1867  with  the  significant  invention  of  dynamite  [Ref.  1],  however  the 
scientific  understanding  of  the  physics  behind  the  detonation  phenomena  was  lacking  and 
required  substantial  research  to  develop  the  technology  to  the  level  existing  today. 
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In  the  early  1870’s,  researchers  such  as  F.  A.  Abel  and  M.  Berthelot  began  to 
relate  the  strength  of  an  explosion  to  how  it  was  initiated  and  its  propagation  velocity.  It 
was  hypothesized  that  some  form  of  mechanical  or  vibrational  shock  waves  must  initiate 
the  detonation  and  that  this  wave  must  be  sustained  throughout  the  event.  M.  Berthelot 
and  P.  Vieille  furthered  this  idea  by  analyzing  the  first  measurements  of  detonation  wave 
velocity  conducted  by  F.  A.  Abel  in  1874.  The  results  were  substantiated  in  similar 
experiments  in  the  early  1880’s.  Bethelot  and  Vieille  concluded,  based  on  numerous 
tests  of  different  fuels  and  oxidizers  at  different  equivalence  ratios,  that  the  detonation 
velocity  is  uniform  and  only  dependent  on  the  fuel  and  its  mixture  ratios.  [Ref.  1] 

E.  Mallard  and  H.  Le  Chatelier  found  another  significant  property  of  the 
detonation  phenomenon  in  1883.  Using  a  new  device  called  a  drum-camera,  they  were 
able  to  demonstrate  that  a  deflagration,  under  the  right  conditions,  will  transition  into  a 
detonation  wave.  This  new  information  led  to  their  experimental  proof  that  the 
detonation  process  can  be  seen  as  a  rapid  adiabatic  process  whose  energy  drives  the 
detonation  wave.  It  was  also  noted  that  the  detonation  wave  velocity  was  on  the  same 
order  as  the  sound  velocity  of  the  products  of  the  combustion.  This  led  to  the  formulation 
that  there  is  an  inherent  relation  between  the  chemistry  of  the  reactants,  the  conditions  in 
which  they  are  ignited,  and  the  detonation  properties  they  exhibit.  [Ref.  1] 

2.  Formulation  of  the  Chapman-Jouget  Theory 
The  early  discoveries  and  ideas  produced  in  the  1880’s  led  to  several  of  the  top 
thermodynamicists  and  chemists  of  the  day  to  become  interested  in  the  new  detonation 
phenomenon.  In  1 890,  V.  A.  Michelson,  working  through  the  Rankine  Theory,  showed 
the  detonation  pressure  as  a  function  of  the  detonation  velocity  and  the  reactant's  heat  of 
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reaction,  and  that  through  the  system  of  equations  there  are  two  possible  solutions.  This 
is  the  first  insight  in  that  reactants,  at  different  conditions,  will  bum  naturally  around  two 
distinct  conditions.  He  also  noted  that  there  exists  a  convergence  of  pressures  at  the 
upper  point,  which  correlated  to  the  detonation  process.  This  was  later  proven  by  the 
work  of  D.  L.  Chapman  and  E.  Jouget.  [Ref.  1] 

Chapman-Jouget  theory  was  recognized  as  the  relationship  between  velocities  of 
combustion  wave  processes  and  the  pressures  at  which  they  occur.  It  also  is  the 
realization  that  there  are  two  regions  at  which  a  combustion  process  may  occur.  The 
theory  is  named  for  two  individuals  who  worked  independently  and  demonstrated  two 
distinct  properties  for  combustion  events.  D.  L.  Chapman  published  his  findings  in  1899 
that  stated  there  exists  a  minimum  velocity  in  which  a  detonation  can  occur  and  it  is 
thermodynamically  tied  to  the  properties  of  the  burned  gas.  The  second  part  of  the  theory 
was  provided  by  the  work  of  E.  Jouget.  He  worked  from  1901-1905  and  established  the 
relation  that  the  detonation  wave  velocity  is  equal  to  the  sound  velocity  of  the  burned  gas 
in  which  it  propagated.  He  verified  this  result  by  comparing  computed  results  with 
experimental  data  found  by  several  of  his  predecessors.  J.  L.  Crussard  validated  the  C-J 
theory  in  1907  by  relating  the  two  specific  combustion  pressure  points  on  the  Hugoniot 
curve  (pressure-specific  volume  adiabat).  The  two  solutions  of  the  C-J  theory  can  be 
seen  on  the  well-publicized  graph  shown  in  Figure  [1-1].  The  upper  C-J  point  is  related 
to  a  detonation  and  the  lower  C-J  point  showing  the  deflagration  region.  [Ref.  1]  The 
relationship  of  the  Chapman-Jouget  points  with  the  Hugoniot  curve  and  a  description  of 
the  distinct  solution  regions  will  be  looked  at  in  detail  in  Chapter  II.  Research  continued 
into  the  twentieth  century  as  Von  Neumann,  Doring  and  Zeldovich  further  modeled  the 
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detonation  wave  as  a  shock  wave  followed  closely  by  the  combustion  process  of  the 
reactants. 


Figure  [1-1].  Chapman- Jouget  Solutions  on  the  Hugoniot  Curve.  [From  Ref.  2] 

3.  Past  Applications  of  Detonations 

With  the  physics  of  the  detonation  wave  explained  using  C-J  theory,  many 
applications  could  be  modeled  and  utilize  the  energy  release  rate  of  detonations.  Before 
the  1890’s,  detonations  were  commonly  used  in  explosions,  fireworks,  and  in  dynamite. 
However,  in  these  applications  the  detonation  was  not  controlled  or  readily  understood. 
The  first  unsuccessful  attempt  to  utilize  the  thermal  efficiency  of  a  detonation  in  a 
propulsion  system  was  seen  in  the  German’s  “Buzz  Bombs”  that  created  havoc  among 
the  Allied  cities  in  World  War  II.  Another  attempt  at  applying  pulsed  combustion  in  a 
propulsion  system  occurred  in  the  mid  twentieth  century  in  the  United  States  under 
“Project  Squid”.  This  program  studied  the  application  of  pulsed  engines  in  both  military 
and  civilian  combustion  engines.  These  pulsejet  engines  were  shown  not  to  be  detonation 
engines  since  their  combustion  processes  occurred  with  a  subsonic  flame  speed  in  a 
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"tuned-port"  combustion  chamber.  These  efforts  were  later  abandoned  since  it  was 
determined  that  the  propulsion  system  did  not  have  a  high  overall  efficiency.  [Ref.  3] 

C.  CURRENT  DETONATION  RESEARCH 

There  is  a  lot  of  promising  research  in  detonation  fundamentals  and  applications 
by  government  labs,  commercial  enterprises,  and  academic  facilities.  The  concept  of 
utilizing  detonations  for  propulsion  was  renewed  in  the  1960’s  with  propulsion  concept 
studies  for  new  high-speed  unmanned  vehicles.  Since  that  time,  there  have  been  many 
studies,  experiments  and  simulations  run  by  several  independent  organizations.  The  most 
documented  and  perhaps  most  well  known  of  these  studies  include  the  research  of  Adroit 
Systems  Incorporated  and  the  Office  of  Naval  Research  funded  academic  consortium. 
Each  of  these  independent  groups  has  aggressive  test  plans  in  hopes  to  gain  insight  in  the 
intricacies  of  effectively  creating  and  controlling  detonations  for  a  useful  propulsion 
system  design. 

Adroit  Systems  Incorporated  (ASI)  has  taken  a  vocal  lead  in  the  commercial 
development  of  detonation  technology.  T.  Bussing  and  G.  Pappas  of  ASI  published  a 
technical  paper  for  the  American  Institute  of  Aeronautics  and  Astronautics  conference  in 
January  1994  that  summarized  the  technical  feasibility  of  designing  and  building  a  Pulse 
Detonation  Engine  (PDE).  This  paper  outlined  some  basic  physics  of  detonations  and 
arguments  supporting  future  research  and  development  of  a  PDE  [Ref.  4],  Since  that 
time  ASI  has  developed  and  patented  many  pulse  detonation  technologies,  such  as  using 
multiple  combustors  with  a  rotary  valve.  They  have  also  invested  in  a  Pulse  Detonation 
Rocket  Engine  (PDRE)  concept  and  using  detonations  for  material  synthesis  and  intense 
sound  level  generation  [Ref.  5]. 
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The  United  States  Department  of  Defense  is  investigating  detonation  technology 
for  use  in  future  high-speed  missile  and  aircraft  propulsion  systems  or  as  a  non-lethal 
weapon.  The  Office  of  Naval  Research  (ONR)  has  organized  an  integrated  research  team 
in  hopes  to  understand  the  requirements  needed  for  PDE  development.  This  team  consists 
of  six  prominent  university  laboratories  and  the  Navy’s  research  facilities  including  the 
Naval  Research  Laboratory  (NRL),  Naval  Postgraduate  School  (NPS)  and  the  Naval  Air 
Warfare  Center  (NAWC).  The  combined  work  of  these  institutions  will  cover  PDE 
research  from  fundamental  detonation  studies  to  multi-cycle  operation  and  cycle  analysis. 
Figure  [1-2]  is  a  slide  outlining  the  detonation  research  roadmap  and  the  responsibilities 


of  the  participating  institutions. 

XNRlS  ONR  PULSE  DETONATION  ENGINE  RESEARCH  ROADMAP 


OtfKKr  Participant* 

—  ssm 

-*  International 


Figure  [1-2].  ONR  Pulse  Detonation  Research  Roadmap.  [From  Ref.  6] 

This  thesis  will  explore  some  of  the  topics  NPS  is  tasked  to  research  in  part  one  of  this 
roadmap.  Specifically,  it  researches  the  ignition  energy/power  requirements  for 
minimization  of  XDdt  as  a  function  of  ignition  location  and  equivalence  ratio  for  an 
ethylene  and  oxygen  fuel  mixture.  [Ref.  6] 
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II.  PHYSICS  OF  A  DETONATION  WAVE 


A.  INTRODUCTION 

In  order  to  begin  a  detailed  analysis  of  detonation  waves,  one  must  first  look  into 
the  subject  of  combustion.  Combustion  is  defined  as  an  exothermic  chemical  reaction 
between  a  fuel  and  an  oxidizer  that  once  initiated  can  sustain  itself  as  long  as  the 
ingredients  in  the  proper  proportions  are  present.  However,  not  all  combustion  events  are 
the  same.  The  velocity  at  which  a  combustion  wave  propagates  through  a  mixture  is  an 
accurate  measure  of  the  strength  or  violence  of  the  event.  This  velocity  is  dependent  on 
several  factors  including  the  mixture  composition,  pressure,  temperature,  and  the 
geometry  of  the  volume  where  the  combustion  occurs.  The  following  sections  describe 
three  distinct  combustion  events  and  how  they  are  initiated,  sustained,  and  terminated. 

B.  DEFINITIONS 

1.  Deflagration 

Deflagration  is  the  combustion  process  most  commonly  utilized  in  current  flight 
propulsion  systems  such  as  ramjets,  turbojets  and  rocket  engines.  Some  practical 
examples  of  this  type  of  combustion  include  a  burning  match,  a  campfire,  or  a  flame 
moving  along  an  open  gasoline  spill.  A  closer  look  reveals  that  a  deflagration  is  a 
combustion  wave  propagating  at  a  subsonic  speed.  The  deflagration  flame  speed  is  a 
function  of  pressure,  temperature,  and  turbulence  of  the  reactants,  and  the  permittable 
range  is  correctly  predicted  by  the  C-J  theory  falling  somewhere  near  the  lower  C-J  point. 
This  speed  is  heavily  dependent  on  the  chemical  makeup,  mass  diffusion  rates  and 
thermal  transfer  rates  of  the  reactants.  Typical  wave  speeds  for  a  deflagration 
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combustion  wave  range  from  l-30m/s.  Though  deflagration  is  the  most  common  and  most 
utilized  combustion  process,  it  is  theoretically  not  the  most  efficient  thermodynamic  path 
for  combustion  to  occur.  This  is  because  the  entropy  of  the  resulting  gases  is  maximized, 
which  reduces  the  amount  of  work  available.  [Ref.  2] 

2.  Explosion 

An  explosion  can  be  a  combustion  or  non-combustion  event  and  occurs  when  a 
reaction  releases  energy  at  a  much  greater  rate  than  the  surrounding  environment  can 
absorb.  This  exothermic  reaction  rate  increases  exponentially  with  the  subsequent 
increase  in  temperature  and  pressure.  This  causes  the  combustion  event  to  go  quickly  out 
of  control.  The  volume  of  the  explosion  will  expand  due  to  the  increase  of  pressure  and 
break  any  boundaries  placed  on  it.  Even  though  the  explosion  is  powerful  and  occurs 
very  fast,  the  combustion  event  itself  occurs  as  a  deflagration  wave  as  it  travels  through 
unbumed  reactants.  This  becomes  an  important  distinction  between  the  explosion  and 
the  detonation  wave.  [Ref.  3] 

3.  Detonation 

A  detonation  is  a  combustion  event  that  happens  when  special  conditions  exist.  It 
can  be  described  as  a  supersonic  combustion  event  that  propagates  at  high  velocities  and 
produces  a  violent  and  rapid  combustion  of  the  reactants  due  to  the  strong  shock  wave 
leading  the  detonation.  The  mechanics  of  the  detonation  will  be  derived  and  discussed  in 
detail  in  Section  C  of  this  chapter. 

A  qualitative  distinction  between  detonations  and  deflagrations  is  a  good  format 
to  view  the  extreme  differences  between  the  two.  The  reference  frame  for  this  analysis 
will  be  a  one-dimensional  stationary  combustion  wave.  The  properties  used  to  describe 
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the  event  will  include  ratios  of,  reactant  velocity  ( u ),  sonic  velocity  (c),  density  (p), 
temperature  (7)  and  pressure  (p).  Subscripts  will  delineate  the  burned  and  unbumed  side 
of  the  wave.  Figure  [2-1]  is  a  schematic  diagram  of  the  reference  frame  used  for  the 
analysis. 


Stationory  Combustion  Wove 

v///////////////////////////^/////////////// 


(Unburned) 
u. 


'I 


p».Ti.pi 


(Burned) 


u„ 


P2,  T2,  P2 


77777Z!7777^P77^77V77777777777?7}777^7777777? 

Figure  [2-1],  Schematic  of  Stationary  1-D  Combustion  Wave.  [From  Ref.  2] 

A  qualitative  comparison  between  a  deflagration  and  detonation  is  given  in  Table  [2-1]. 

Note  that  the  detonation  pressures,  velocities  and  density  ratios  are  significantly  larger 

than  the  deflagration  process.  These  results  depict  dramatic  differences  between  the  two 

combustion  events  and  lead  to  questions  as  to  how  it  happens,  and  what  can  be  gained 

from  it. 


Table  [2-1],  Qualitative  Difference  Between  Detonation  and  Deflagration.  [From  Ref.  21 


Properties 

Deflagration 

Detonation 

U\/C\ 

0.0001-0.03 

5-10 

Uilu\ 

4-6  (acceleration) 

0.4  -  0.7  (deceleration) 

Pllpx 

=  0.98  (slight  expansion) 

13  -55  (compression) 

T2/T] 

4-16  (heat  addition) 

8-21  (heat  addition) 

P2/P1 

0.06  -  0.25 

1.7 -2.6 

C.  DETONATION  WAVE  THEORY 

1.  Thermodynamics  of  a  Detonation  Wave 

Any  detailed  study  into  the  origins  and  properties  of  a  detonation  wave  must 
begin  with  a  derivation  and  analysis  of  the  Hugoniot  curve.  This  relationship  between 
enthalpy  (h),  pressure,  and  density  of  the  gases  in  a  combustion  event  is  directly  related 
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to  the  various  combustion  conditions.  The  Chapman-Jouget  points  discussed  in 
Chapter  1  are  the  two  physical  solutions  to  the  Hugoniot  relation  for  the  constant  area 
geometry  given,  and  will  the  basis  for  analyzing  the  detonation  phenomena. 

The  governing  equations  of  a  thermodynamic  process  must  be  derived  from  the 
basic  conservation  equations  of  physics.  These  equations  are: 

Continuity  equation:  ~*dx^ =  ^  [2-1] 


Conservation  of  momentum: 


du  d(pu )  dp 
pu — +u  = — — 

dx  dx  dx 


Conservation  of  energy: 


pu 


d_ 

dx 


(  ..2\ 


h  + 


u 


=  -~r  {qCOnd) 
dx 


Where  the  enthalpy  ( h )  and  the  heat  added  {qCOnd)  to  the  system  is  defined  by 


[2-2] 

[2-3] 


h  =  CpT  +  h° 

[2-4] 

q  =  h;-h; 

[2-5] 

ldT 

cond  ^  1 

dx 

[2-6] 

This  derivation  assumes  steady  one-dimensional  flow,  with  no  external  heat  added  or 
rejected,  negligible  interdiffusion  effects,  no  viscous  effects,  and  occurs  within  a  constant 
area  combustor.  By  analyzing  the  detonation  in  this  form,  it  can  be  viewed  as  a 
supersonic  shock  wave  with  calculable  properties  in  front  of  and  behind  the  wave.  By 
integrating  equation  [2-1]  the  concept  of  a  constant  mass  flow  rate  (m )  is  revealed  by 


jd(pu)  =  j Odx : 


pu  =  const.  =  m 


[2-7] 


Substituting  equation  [2-1]  into  the  momentum  equation  and  simplifying,  equation  [2-2] 
becomes 
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[2-8] 


Au,2 

dx 


[pu2+p}  =  0 


Integrating  equation  [2-8]  gives  an  alternate  momentum  equation  to  work  with. 

pu2  +  p  =  const.' 

By  similar  analysis,  the  energy  equation  [2-3]  becomes 


[2-9] 


( 


u 


+  h°  +  ^ 


dT 

-A —  =  const." 
dx 


[2-10] 


Referring  back  to  Figure  [2-1]  and  knowing  that  the  change  in  temperature  with  respect 
to  position  ( dT/dx )  in  front  of  and  behind  the  detonation  wave  is  equal  to  zero,  specific 
conservation  equations  can  be  derived.  These  equations  relate  the  thermodynamic  and 
fluid  dynamic  properties  between  the  two  regions  of  the  reference  frame.  These 
equations  are: 


P\u\  -  Piui  =  const.  =  m 

[2-11] 

Pi  +P\U\  =  P2+P2Ul 

[2-12] 

2  2 

CpT]+-f+q  =  CpT2+^- 

[2-13] 

Using  equation  [2-4]  and  [2-5],  equation  [2-12]  can  be  written  as 

u 2  u\ 

=  [2-14] 


The  fourth  and  final  equation  needed  for  deriving  the  Hugoniot  relation  is  found  in  the 

assumption  that  the  gases  in  both  the  burned  and  unbumed  regions  behave  like  a  perfect 

p  =  pRT  [2-15] 

gas.  The  perfect  gas  law  is  defined  for  both  regions  and  R  is  the  specific  gas  constant  for 

the  reactants.  Equations  [2-1 1],  [2-12],  [2-14]  and  [2-15]  are  the  basis  set  of  equations 
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used  to  solve  the  five  unknowns  of  the  system,  uu  u2,  p2,  T2,  and  p2.  One  equation  with 
two  unknowns  (p2,/?2)  can  be  formed  by  manipulating  the  four  basis  equations  in  the 
following  manner: 


Pi  P\  ~  P\u\  P1P1 


(PlUlf  (p2Ul)2 
Pi  Pi 


Substituting  equation  [2-7]  in  to  [2-16]  yields 


_  {PiMi ) 

Pi  Pi 


yPi 


Pi) 


Therefore, 


2  Pl~  P\  *2 
Piui  =ii  =m 


Pi  Pi 


[2-16] 


[2-17] 


[2-18] 


Equation  [2-18]  is  the  Rayleigh-line  relation,  which  quantitatively  describes  the  heat 
addition  process  in  relation  to  velocity  (w),  pressure  (p),  and  density  (p)  in  a  gas. 
Rearranging  the  Rayleigh-line  relation  in  terms  of  Mach  number  (M)  is  a  more  useful 
way  to  describe  the  characteristics  of  the  thermodynamics.  This  is  done  with  the 
following  approach: 


a  |  o 
ll 

[2-19] 

Where,  c  =  ^RT 

[2-20] 

Substituting  equation  [2-15]  gives, 
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Manipulating  equation  [2-18]  and  substituting  in  equation  [2-21],  the  Rayleigh-line 
relation  becomes: 


Y_  Plu\ 

r  PiPi 


^-1 

=  - =  jM? 

1-A 

Pi 


[2-22] 


Knowing  that, 


Cp  ~CV  =  R=>  Cp  =  • 


R 


(_Y_' 

<r- 1. 


[2-23] 


/ 


equation  [2-13]  can  be  rewritten  as, 


Y 


( 


r- 1 


p±_p± 

/^2  Z7!  y 


~{u2x~ul)=q 


[2-24] 


Using  equation  [2-12]  to  create  an  expressions  for  u\  and  u\,  and  substituting  into 
equation  [2-24],  the  relation  becomes, 


r 


r- 1 


Pi__P\ 
Pi  P\  j 


— — —  +  —u\  +2L-Pi 


=  q 


[2-25] 


\  P\  P\  P 2  Pi  J 

Combining  equations  [2-1 1]  and  [2-25]  gives  the  final  relation,  equation  [2-26],  between 
heat  addition  and  the  gases  initial  and  final  pressures  and  densities. 


Y 


r- 1 


\Pi  P\ 


i/  {i  n 

—z{P2-pA  —  +  —  : 
2  vZ7!  Pi  J 


[2-26] 


This  relation  can  also  be  written,  with  some  manipulation  of  equation  [2-4],  [2-5]  and 
[2-15]  and  substituting  into  equation  [2-26]. 


hi~K  =\{pi-pi—+— 

2  Iz7,  Z72 . 


[2-27] 
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These  final  two  equations  are  forms  of  the  Rankine-relation,  but  formally  named  for 
Hugoniot  who  derived  them  and  fit  the  various  combustion  conditions  to  the  pressure 
versus  specific  volume  ( 1 / p)  curve. 

The  Hugoniot  curve  is  a  plot  that  describes  the  different  possible  combustion 
conditions.  These  combustion  conditions  include  various  strengths  of  deflagrations  and 
detonations  dependent  upon  the  pressure  and  specific  volume  conditions  at  which  the 
event  is  occurring.  The  Chapman-Jouget  points  described  earlier  can  be  plotted  on  this 
graph  and  prove  to  be  the  boundaries  for  strong  and  weak  combustion  events. 

Figure  [2-2]  is  the  Hugoniot  curve  along  with  a  depiction  of  the  different  sections  and 
boundaries  of  combustion  events.  The  origin  of  the  Hugoniot  curve  is  labeled  A  and  is 
constructed  by 


Figure  [2-2].  Hugoniot  Curve  Depicting  Different  Combustion  Regions.  [From  Ref.  2] 

intersecting  the  tangents  to  the  upper  ( U)  and  lower  ( L )  C-J  points.  The  combustion 

regions  for  both  detonations  and  deflagrations  can  be  described  as  both  strong  and  weak, 

as  depicted  by  regions  I  through  IV  in  Figure  [2-2].  Region  V  is  a  portion  of  the  curve 
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that  is  physically  impossible  in  that  it  needs  uj  to  be  imaginary  to  satisfy  the  Rayleigh- 
line  relation.  However,  the  region  of  interest  for  studying  detonation  transitions  lies  near 
the  upper  C-J  point.  It  can  be  shown,  and  will  be  seen  in  the  discussion  of  the  detonation 
wave  structure,  that  a  detonation  tends  to  naturally  occur  at  this  critical  point  on  the 
Hugoniot  curve.  [Ref.  2] 

2.  Structure  of  a  Detonation  Wave 

A  detonation  wave  is  sustained  through  a  complex  balance  between  a  shock  wave 
and  the  combustion  event  behind  it.  This  symbiotic  relationship  is  best  understood 
through  a  one-dimensional  analysis  of  a  detonation  wave.  Three  scientists  completed  the 
analysis  independently  in  the  early  1940's.  Zel'dovich,  von  Neumann,  and  Doring 
believed  the  detonation  wave  could  be  viewed  as  three  distinct  regions  whose  widths  are 
heavily  dependent  on  the  equivalence  ratios  and  kinetics  of  the  gaseous  mixture.  Figure 
[2-3]  is  a  depiction  of  the  thermodynamic  properties  in  the  three  regions  of  the  commonly 
named  ZND  detonation  wave  structure.  The  first  region,  the  shock  wave,  has  a  width  of 
just  a  few  Angstroms,  yet  delivers  a  tremendous  amount  of  energy  into  the  unbumed 
reactants  seen  as  increases  in  pressure,  density,  and  temperature.  The  dramatic  increase 
in  the  thermodynamic  properties  increases  the  chemical  reaction  rates  and  accelerates  the 
energy  release  phase  of  the  wave  structure.  The  deflagration  region  is  made  of  two  zones 
that  describe  the  thermodynamics  of  combusting  the  reactants.  The  first,  which  is  known 
as  the  Induction  zone,  has  a  relatively  short  width  in  which  the  chemical  reaction  is 
beginning  but  not  yet  impacting  the  thermodynamic  properties.  This  shifts  to  the 
Reaction  zone  when  the  reaction  rate  increases  exponentially,  driving  temperatures  up 
and  stabilizing  pressure  and  density  to  their  final  equilibrium  values.  These  three  zones 
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Shock 


Figure  [2-3].  Physical  Properties  of  the  1-D  Detonation  Wave  Structure.  [From  Ref.  2] 
all  occur  within  a  distance  on  the  order  of  a  centimeter  and  each  is  dependent  on  the  next 
to  sustain  the  detonation  wave. 

D.  THERMODYNAMIC  EFFICIENCIES 

1.  Increased  Thermal  Efficiency  of  a  Detonation  Wave 
The  motivation  for  utilizing  detonations  as  a  means  of  propulsion  is  the  possibility 
of  increased  efficiency.  This  is  not  surprising  because  nearly  all  significant 
breakthroughs  in  combustion  engines  have  come  in  the  form  of  decreasing  the 
thermodynamic  losses  within  the  engine.  Thermodynamically,  the  efficiency  of  any 
system  is  directly  related  to  its  ability  to  minimize  the  entropy  rise  in  the  working  fluid. 
Therefore,  a  careful  analysis  of  the  entropy  budget  must  be  completed  to  justify  further 
research  in  detonations  as  a  means  of  propulsion. 

An  expression  for  the  entropy  (s)  of  a  combustion  event  can  be  derived  from  the 
Hugoniot  relationship  discussed  earlier.  In  order  to  prove  the  increased  efficiency  of  a 
detonation,  the  entropy  at  the  upper  C-J  point  on  the  Hugoniot  curve  must  be  analyzed. 
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Kuo,  in  Principles  of  Combustion,  describes  this  derivation  in  the  following  steps  [Ref. 

2].  To  begin,  the  energy  of  the  system  must  be  defined  and  well  understood.  The 
enthalpy  (h)  is  defined  as 

h  =  e+—  [2-28] 

P 

where  e  is  the  total  internal  energy  within  the  system.  Applying  this  two  the  1-D 
stationary  detonation  wave,  equation  [2-28]  becomes 


*2 -A,  )+(—-—  1  [2-29] 

\Pi  P\  J 

Substituting  equation  [2-29]  into  the  Hugoniot  relation,  equation  [2-27],  and  sorting 
terms  gives 

1  i  i  > 

(e2  ~£\)  —  —  (p 2  —  P\ )+  [2-30] 

2  \Pi  PiJ 

Differentiating, 

1  f  1  1  ^ 

de2  =-(dP2\ - -  d~  (p2  +P\)  [2-31] 

2  V  A  Pi  J  ^  V  P) 

Combining  the  first  and  second  laws,  the  following  expression  can  be  generated  which 
introduces  entropy  into  the  analysis. 

(  1  ^ 

T2ds2  =  de2  +  p2d  —  [2-32] 

\Pt  J 


Substituting  equation  [2-31]  into  [2-32]  and  rearranging  terms,  gives  an  expression  in 
terms  of  a  change  in  entropy  over  a  change  in  specific  volume 
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Pi- Pi 


T  ds2  =  1  [  J _ 1_]  _JPi__PiZ£i_ 

2  d(\/p2)  2{pl  p2J  d(\/p2)  J _ 1_ 

v  Pi  P\ 


[2-33] 


Analyzing  the  Hugoniot  curve  at  the  upper  and  lower  Chapman- Jouget  points,  the  slope 
of  the  tangents  to  the  curve  at  these  points  is 


dp  i 

d{yp2) 


p2 -Pi 
_i _ i_ 

pi  Pi 


[2-34] 


Therefore,  when  equation  [2-34]  is  applied  to  equation  [2-33]  representing  the  conditions 
at  the  upper  and  lower  C-J  points,  it  becomes 


d(VP2) 


[2-35] 


This  result  indicates  that  the  entropy  value  on  the  Hugoniot  curve  at  the  C-J  points  is 
either  a  maximum  or  a  minimum.  In  order  to  determine  what  the  condition  is,  the  second 
derivative  in  entropy  must  be  taken. 


t  \  |  H  \ 1  I  dlp 2 _ 1  dp >2  P2  ~Pi  n 

’  4/a)2L~2U  a)  4l/pJ  f  1  1  )d(l/p2)  f  !  !  Y  [2'36] 


_  1 f  1  1  ]  d2p2 


Pi  P\ 


Pi  Pi 


Isolating  the  second  derivative  in  entropy  and  using  the  condition  outlined  in  equation 
[2-34],  equation  [2-36]  reduces  to 


1(1  1 


diVPif)^  2  T\Pl  PtXdil/pj)^ 


[2-37] 
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A  separate  analysis  of  the  asymptotes  of  the  Hugoniot  curve  shows  the  second  derivative 
of  pressure  in  equation  [2-37]  is  greater  than  zero. 


d2P* 

d{l/p2)2 


>0 


[2-38] 


Using  this  fact  and  knowing  the  difference  between  the  specific  volumes  at  the  upper  and 
lower  C-J  points,  it  can  be  determined  which  point  maximizes  entropy  and  which 
minimizes  it.  For  deflagrations  the  specific  volume  (1/pi)  behind  the  shock  front  is  less 
than  that  behind  it  (I/P2).  Therefore  for  1/pi  <  l/p2 


LowerC-J 


<0 


[2-39] 


which  corresponds  to  a  post  combustion  condition  with  the  highest  entropy  rise  for  a 
specified  amount  of  heat  addition,  at  the  lower  C-J  point.  Conversely,  at  the  upper  C-J 
point,  analogous  to  the  detonation,  the  opposite  holds  true.  At  this  point  1/pi  >  l/p2, 
resulting  in 


UpperC-J 


>0 


[2-40] 


which  corresponds  to  a  post  combustion  condition  with  the  lowest  entropy  rise  for  a 
specified  amount  of  heat  addition.  Figure  [2-4]  is  a  graphical  representation  of  the 
entropy  variation  over  the  entire  Hugoniot  curve.  This  analysis  demonstrates  the 
increased  amount  of  work  available  and  efficiency  of  the  detonation  process  over 
deflagration  combustion  events 
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Figure  [2-4].  Entropy  Variation  Along  the  Hugoniot  Curve.  [From  Ref.  2] 

2.  Engine  Cycle  Analysis  Comparison 
a.  Comparison  Approach 

A  practical  approach  to  demonstrating  the  increased  thermodynamic 
efficiency  of  a  Pulse  Detonation  Engine  (PDE)  is  to  compare  the  two  idealized 
combustion  processes  that  could  be  utilized  for  heat  addition.  This  section  will  perform  a 
cycle  analysis  on  a  constant  pressure  (Cp)  combustion,  currently  used  by  turbojets, 
ramjets,  etc.,  and  a  detonation  (C-J)  cycle,  which  exists  for  a  PDE  system.  Both  cycles 
will  be  analyzed  as  ideal  with  the  understanding  they  will  not  completely  model  the 
performance  limit  of  a  real  world  engine.  The  approach  will,  however,  bound  the  real 
world  benefits  of  a  PDE  cycle.  This  section  will  review  the  goals,  assumptions,  and 
equations  needed  to  outline  the  analysis  for  each  cycle. 

The  ultimate  goal  of  this  analysis  is  to  compare  two  idealized  combustion 
cycles  using  the  First  Law  of  Thermodynamics  and  reveal  the  increased  thermodynamic 
efficiency  of  a  detonation  combustion  cycle.  In  doing  this,  a  baseline  understanding  of 
the  fundamental  differences  between  Cp  and  C-J  cycles  will  be  presented.  This  analysis 
will  also  build  a  solid  thermodynamic  approach  to  each  cycle.  The  result  will  be  a  simple 

thermodynamic  analysis  for  the  Pulse-Detonation  Engine  (PDE)  cycle. 
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The  analysis  involves  several  assumptions  of  the  combustion  parameters 
and  cycle  behavior.  The  working  fluid  that  will  be  used  is  an  Ethylene  (C2H4)  and  air  gas 
mixture  at  stoichiometric  conditions.  The  stoichiometric  condition  is  defined  as  the 
mixture's  fuel  to  air  ratio  that  facilitates  complete  combustion,  leaving  no  fuel  or  oxidizer 
in  the  combustion  products.  A  mixture's  equivalence  ratio  is  based  on  the  reactant's 
actual  fuel  to  air  ratio  divided  by  that  of  the  stoichiometric  condition.  The  mixture  is 
assumed  to  behave  in  accordance  with  the  perfect  gas  law, 

pv  =  RT  [2-41] 

The  mixture  will  be  treated  as  thermally  and  calorically  perfect  with  constant  specific 
heats  through  the  isentropic  portions  of  the  cycle.  This  allows  the  enthalpy  (h)  and 
internal  energy  ( u  *)  to  be  functions  of  temperature  (T)  only, 

h  =  CpjdT  [2-42] 

The  energy  released  by  the  combustion  process  will  be  assumed  to  be  simple  heat 
addition  (qi.2)  based  on  the  type  of  combustion  process  for  the  cycle.  Further 
assumptions  must  be  placed  on  the  cycles  in  order  to  standardize  them  for  a  coherent 
analysis.  Both  cycles  are  started  with  an  initial  condition  equivalent  to  that  of  an 
operating  system  at  an  altitude  of  50,000  ft.,  and  will  be  labeled  as  state  (0).  Next,  both 
the  compression  (0-1)  and  expansion  (2-4)  paths  will  be  modeled  as  isentropic.  The  heat 
addition  path  occurs  between  state  (1)  and  (2)  and  will  be  supplied  by  the  combustion 
process.  The  cycle  will  be  completed  by  a  constant  pressure  cooling  which  is  physically 
seen  as  a  new  load  of  reactants  at  the  initial  condition  along  the  (4-0)  path.  Figure  [2-5] 
is  a  generalized  plot  of  a  combustion  cycle  to  outline  the  four  paths  within  the  cycle. 
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Figure  [2-5].  P-v  Diagram  of  Generalized  Combustion  Cycle. 

It  will  be  shown  by  careful  analysis  that  this  plot  shifts  with  different  combustion  events 
that  correspond  to  differences  in  the  total  work  output  and  efficiency  of  the  cycle.  Table 
[2-2]  outlines  the  state  (1)  parameters  that  will  be  used  for  each  cycle  [Ref.  7],  Each 
cycle  will  assume  an  initial  compression  ratio  of  10:1.  Parameters  at  each  state  of  the 
cycles  will  be  given  during  the  analysis.  These  values  will  make  it  possible  to  calculate 
the  work  and  efficiency  of  that  cycle  in  order  to  adequately  compare  them. 


Table  [2-2],  Standard  Parameters  for  Combustion  Cycle  Analysis.  [From  Ref.  7] 


•  Parameter 

Symbol 

Value 

Units 

y 

1.35 

None 

Specific  Heat 
(constant  pressure) 

cp 

1.04 

kJ/kg 

Standard 

Temperature 

T 

K 

Standard  Pressure 

P 

atm 

Specific  Volume 

V 

mJ/kg 

Entropy 

s 

6.48 

kJ/kgK  | 

The  equations  utilized  for  the  analysis  are  derived  from  the  basic  laws  of 
physics  and  thermodynamics.  As  previously  outlined,  the  conservation  of  mass, 
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momentum  and  energy  along  with  the  second  law  of  thermodynamics  are  required  to 
perform  the  cycle  analysis.  These  equations  in  cycle  notation  are  presented  again  for 
convenience: 


Conservation  of  mass: 

pxux  =  p2u2  =  m 

[2-43] 

Conservation  of  momentum: 

plux+p]  =  p2u\+p2 

[2-44] 

Conservation  of  energy: 

,  u\  .  ul 

K  +-^  +  (h->2  =h2  +  — 

[2-45] 

Or 

r  W]  f  tl2 

U\+P\Vl+Y  +  (h->2=U2+P2V2+Y 

[2-46] 

substituting  in  the  relation, 

h  =  u'  +  pv 

[2-47] 

The  second  law  of  thermodynamics  combined  with  the  Gibb's  free  energy  function  will 
also  be  utilized  to  determine  the  entropy  between  the  phases  of  the  cycles.  It  follows  the 
form, 

s  =  C,  ln(r/7-Rt/ )  -  R 1 n(p/pfc/  )  [2-48] 

Another  relation  that  defines  how  a  combustible  gas  will  behave  in  a  given  cycle  is  that  it 
must  follow  the  Hugoniot  curve  and  Rayleigh  line  relation  defined  earlier  in  Section  C  of 
this  chapter.  The  curves  define  all  the  mathematically  possible  end  states  for  a  gas  with 
the  specified  amount  of  heat  release  and  initial  conditions,  although  only  certain  regions 
of  the  curve  are  attainable.  The  curves  are  hyperbolic  in  form  and,  by  manipulating  the 
three  conservation  equations,  the  velocity  of  the  gas  dependence  can  be  removed.  All  of 
these  relations  must  be  satisfied  for  a  given  combustion  process  to  occur  in  nature,  and 
using  them  will  produce  the  cycle  parameters  required  for  evaluating  cycle  efficiency. 
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Combustion  cycle  efficiency  will  be  calculated  thermodynamically  using 
the  standard  definition, 

W 

Pth  =  [2-49] 

P  addition 

However,  the  total  work  of  the  cycle  must  be  calculated  first.  This  will  be  done  utilizing 
the  conservation  of  energy  principle. 

jSV  =  c \Sq  [2-50] 

By  solving  for  the  total  work  produced  by  the  combustion  cycle,  the  thermal  efficiency 
can  be  solved  by  knowing  the  amount  of  heat  released  from  the  reactants  during  the 
combustion  event.  The  thermal  efficiency  of  both  cycles  will  be  conducted  to  determine 
the  possible  thermodynamic  gains  a  detonation  based  cycle  could  give  in  a  propulsion 
system. 

b.  Deflagration  at  Constant  Pressure  Analysis 

This  analysis  will  step  through  the  procedure  used  to  determine  the  state 
values  through  a  constant  pressure  combustion  cycle.  The  initial  conditions  used  for  all 
calculations  is  outlined  in  Table  [2-2]  at  the  pre-combustion  state  (1).  To  solve  for  the 
initial  state  (0),  the  isentropic  relationship,  shown  by  equation  [2-51],  coupled  with  the 
10: 1  compression  ratio  is  used.  The  temperature  at  state  (0)  is  then  solved  for  utilizing 
the  perfect  gas  law  shown  by  equation  [2-41]. 

Pvr  =  Const.  [2-51] 

To  solve  for  the  post-combustion  state  (2),  the  thermodynamic  code  TEP  was  run  for 
constant  pressure  conditions  with  the  defined  reactant  ratios  [Ref.  8].  The  output  from 
TEP  provided  the  post-combustion  temperature,  y,  Cp,  and  entropy  of  the  mixture. 
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Assuming  the  products  behave  as  a  perfect  gas,  the  specific  volume  for  the  state  can  be 
solved  using  the  perfect  gas  law.  The  properties  at  state  (4),  are  arrived  at  through  an 
isentropic  expansion  and  are  solved  for  using  the  same  isentropic  analysis  used  for  the 
compression  process.  Table  [2-3]  depicts  the  results  of  this  analysis.  A  pressure-specific 
volume  (P-v)  diagram,  shown  in  Figure  [2-6],  displays  the  constant  pressure  combustion 
cycle.  This  diagram  can  be  used  to  compare  the  constant  pressure  cycle  with  the  C-J 
detonation  cycle. 


Table  [2-3].  State  Parameters  for  a  Constant  Pressure  Combustion  Cycle. 


State 

Pressure  (atm) 

Temperature 

(K) 

v  (m3/Kg) 

0 

0.1 

164.9 

5.3 

6.5 

1 

1.0 

300.0 

0.8 

6.5 

2 

1.0 

2353.0 

6.2 

9.6 

4 

0.1 

1361.1 

39.3 

9.6 

Figure  [2-6].  P-v  Diagram  of  Constant  Pressure  Combustion  Cycle. 


The  constant  pressure  combustion  cycle's  thermal  efficiency  calculation 

requires  the  heat  added  by  the  combustion  to  be  determined.  The  difference  between  the 

pre-combustion  and  post-combustion  enthalpies  were  calculated  to  determine  this  value. 
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The  conservation  of  energy  equation  shows  that  this  difference  is  equal  to  the  heat  added 
to  the  reactants  from  combustion.  Equation  [2-52]  depicts  this  calculation  and  is  derived 
by  manipulating  equations  [2-42]  and  [2-45]  and  assuming  no  change  in  the  reactants 
velocity. 

0>2^2  -CnTt  =  hj  =  2  [2-52] 

From  this  calculation,  the  heat  added  by  combustion  (qi2)  is  equal  to  3105.1  (KJ/Kg-K). 
The  heat  released  to  the  atmosphere  is  calculated,  assuming  a  constant  post-combustion 
Cp,  using  equation  [2-53]. 

Q>2^2  —  =  ?4->0  [2-53] 

The  net  work  (Wnet)  of  the  cycle  is  then  calculated  to  be  1749.2  (KJ/Kg-K).  Using 
equation  [2-49],  the  thermal  efficiency  (pth)  of  the  idealized  constant  pressure 
combustion  cycle  is  calculated  to  be  41.9%. 

c.  Detonation  Analysis 

This  analysis  will  step  through  the  procedure  used  to  determine  the  state 

values  through  a  C-J  combustion  cycle.  In  order  to  adequately  compare  the  C-J  cycle  to 

the  constant  pressure  cycle,  the  same  initial  pre-combustion  compression  ratios  were 

used.  The  properties  for  states  (1)  and  (2)  were  found  utilizing  the  same  procedure 

outlined  in  Section  D.2.c.  To  solve  for  the  post-combustion  state  (3),  the  thermodynamic 

code  TEP  was  run  under  detonation  conditions  with  the  defined  reactant  ratios  [Ref.  8]. 

Again,  assuming  the  products  behave  as  a  perfect  gas,  the  specific  volume  for  the  state 

was  solved  for  using  the  perfect  gas  law.  The  properties  at  state  (4),  were  determined 

through  an  isentropic  expansion  and  are  solved  for  using  the  same  isentropic  analysis 

used  for  the  compression  process.  Table  [2-4]  depicts  the  results  of  this  analysis.  A 
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pressure-specific  volume  (P-v)  diagram,  shown  in  Figure  [2-7],  displays  the  C-J 
detonation  cycle  and  can  be  used  to  compare  the  cycle  with  the  Cp  combustion  cycle. 


Table  [2-4],  State  Parameters  for  a  C-J  Detonation  Cycle. 

I  I  T DTnt\orofnrD  I  1 1  />\  I 


State 

Pressure  (atm) 

Temperature 

(K) 

v  (mJ/Kg) 

entropy  (5) 
(KJ/Kg-K) 

0 

0.1 

164.9 

5.3 

6.5 

1 

1.0 

300.0 

0.8 

6.5 

2 

18.2 

2926.0 

0.4 

9.2 

4 

0.1 

1022.9 

27.2 

~  9.2 

The  detonation  cycle's  thermal  efficiency  was  determined  by  analyzing  the 
differences  in  the  enthalpies  in  the  pre  and  post-combustion  states.  The  thermodynamic 
conservation  of  energy  shows  that  this  difference  is  also  equal  to  the  heat  added  to  the 
working  fluid  from  the  combustion  process.  Equation  [2-54]  is  derived  by  manipulating 
equations  [2-42]  and  [2-45].  However,  in  a  detonation  the  velocities  of  the  reactants 
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contain  a  substantial  amount  of  energy  and  cannot  be  ignored  for  the  calculation.  For  an 
assumed  stationary  combustion  wave,  the  equation  becomes: 

(y  -r)2  V2 

CP2T2  A ■*YL-CnTl -~f  =  h}  =  qM  [2-54] 

This  calculation  determines  the  heat  added  by  combustion  (q^)  to  be  equal  to  2960.7 
(KJ/Kg-K).  The  heat  released  to  the  atmosphere  was  calculated  assuming  a  constant 
post-combustion  Cp  and  using  equation  [2-55]. 

CP1Ti  -  CPiT0  =  qA_^  [2-55] 

With  these  values  known,  the  net  work  (Wnet)  of  the  cycle  was  then  calculated  to  be 
1758.2  (KJ/Kg-K).  By  equation  [2-49],  the  thermal  efficiency  (t|th)  of  the  idealized 
constant  pressure  combustion  cycle  equals  65.6%. 

3.  Thermodynamic  Analysis  Results 

The  analysis  carried  out  in  the  previous  two  sections  analyzed  two  separate 
combustion  processes  at  the  same  initial  conditions.  The  results  of  the  analysis  depict  a 
substantial  difference  in  the  overall  thermal  efficiency  of  the  two  cycles.  Table  [2-5]  is  a 
summary  of  the  important  efficiency  parameters,  t|th,  entropy  change  (As),  and  the 
percent  difference  between  the  two  cycles. 


A  graphical  display  of  the  difference  in  the  net  work  can  be  seen  in  the  integration 


under  the  P-v  diagrams  for  both  cycles.  The  C-J  combustion  process  gains  its  advantage 


in  the  nearly  constant  volume  combustion  that  drives  the  pressure  18  times  that  of  the 
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constant  pressure  combustion.  Figure  [2-8]  displays  these  differences,  which  can  be  seen 
by  comparing  the  areas  under  each  curve,  which  represents  the  net  work.  Another,  more 
classical,  way  to  compare  the  two  cycles  is  to  look  at  the  net  efficiency  gain  for  each. 
Following  the  theory  outlined  in  Section  D.l  of  this  chapter,  the  C-J  detonation  cycle 
results  in  a  significantly  lower  entropy  rise,  which  results  in  more  work  available.  For 
this  reason  the  Pulse  Detonation  Engine  concept  is  a  theoretically  attractive  alternative 
for  a  high-speed  long-range  propulsion  system. 


Figure  [2-8].  P-v  Diagram  Comparison  of  C-J  Detonation  vs.  Cp  Combustion  Cycle. 


E.  PHYSICAL  MECHANICS  OF  A  DETONATION 
1.  Deflagration-to-Detonation  Transition 

The  goal  of  this  research  is  to  determine  conditions  for  minimum  deflagration-to- 
detonation  transition  (Xddt)  distance  in  a  simple  tube  combustor  with  one  open  end.  In 
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order  to  quantitatively  measure  this  parameter,  the  concept  of  transitioning  to  a 
detonation  wave  must  be  understood.  Kuo  describes  the  sequence  of  events  leading  to  a 
fully  developed  detonation  wave  in  four  generalized  steps. 

•  First,  laminar  flames  of  an  enclosed  deflagration  produce  compression  waves 
that  coalesce  to  form  a  shock  front  in  the  unbumed  reactants. 

•  This  shock  front  creates  a  turbulent  reaction  zone  within  the  flame  front, 
which  eventually  causes  an  "explosion  within  an  explosion"  behind  the  shock 
front. 

•  This  "explosion"  produces  strong  shock  waves  in  both  directions  and 
oscillations  between  them,  called  transverse  waves.  The  rearward 
propagating  shock  is  called  a  retonation  while  the  forward  shock  is  known  as 
a  superdetonation. 

•  Interactions  between  all  of  these  events  coupled  with  the  reaction  zone  of  the 
primary  flame  front  combine  to  create  a  steady  state  supersonic  C-J 
detonation  wave. 

Figures  [2-9]  and  [2-10]  are  depictions  of  this  developing  detonation  wave  process 
outlined  above.  [Ref.  2] 


Spherical 

Shock 


Figure  [2-9].  "Explosion  Within  Explosion"  in  a  Developing  Detonation.  [From  Ref.  2] 
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Retonotion  Compression 

Wave  Wave 


Figure  [2-10],  Components  of  a  Transitioning  Detonation  Wave.  [From  Ref.  2] 
Experiments  measuring  the  wall  static  pressure  of  the  combustion  wave  front  as  it 
propagates  down  the  combustor  tube  show  that  it  is  possible  to  determine  where  the  wave 
transitions  from  a  deflagration  to  a  detonation.  Figure  [2-1 1]  is  a  depiction  of  a  typical 
transitioning  detonation  wave.  Note  the  gradual  rise  in  pressures  and  shock  front  speed 
in  the  first  two  transducers.  This  is  indicative  of  the  coalescing  of  shock  waves  in  the 
transitioning  process.  Referring  back  to  the  ZND  wave  structure  of  a  detonation,  Figure 
[2-5],  the  fully  developed  detonation  wave  will  show  a  sharp  and  drastic  increase  in 
pressure  as  it  passes  the  transducer.  This  can  be  seen  in  the  wave  trace  of  transducer  P5 
shown  in  Figure  [2-11].  That  is  the  position  at  which  the  wave  is  determined  to  be  fully 
transitioned.  The  retonation  phenomena  can  be  seen  in  the  subsequent  pressure  peaks  of 
P6  and  P7  as  the  retonation  propagates  back  over  the  transducers  and  reflects  off  the  head 
wall  of  the  combustor. 
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Figure  [2-1 1].  Pressure  Trace  of  a  Transitioning  Detonation  Wave. 


By  analyzing  the  shapes  and  speeds  of  the  wave  at  different  points  along  the  combustor, 
it  is  clear  that  the  wave  transitioned  into  a  detonation  between  P4  and  P5  transducers. 
Once  transitioned,  the  wave  continues  to  propagate  at  the  steady-state  C-J  detonation 
velocity. 


2.  Detonation  Wave  Velocity 

The  key  indication  that  a  combustion  event  is  a  detonation  is  to  measure  its  wave 

velocity  and  compare  it  with  theoretical  values.  Detonation  wave  velocities  can  reach 

velocities  of  1500-3000  meters  per  second  dependent  upon  the  chemical  properties  of  the 

reactants.  This  detonation  velocity,  once  reached,  generally  remains  constant  throughout 

the  combustor.  This  fact  is  predicted  by  Chapman-Jouget  theory  of  detonations  and 

proven  in  well-documented  experiments.  This  was  also  seen  in  the  experiments 

conducted  for  this  research.  The  detonation  velocity  (Vdet)  was  determined  using  the 

32 


static  pressure  traces  of  the  wave  as  it  traveled  through  the  combustor.  The  rate  equation 
was  used  for  this  calculation  by  timing  the  wave  as  it  passed  over  a  known  distance 
between  transducers.  Figure  [2-12]  depicts  a  detonation  wave  trace  as  it  propagates  over 
pressure  transducers  spaced  2  inches  apart.  This  fact  along  with  a  measured  time 
between  them,  allows  the  experimental  Vdet  to  be  calculated  to  be  2822±180  meters  per 
second  for  this  particular  run. 


Figure  [2-12].  Detonation  Wave  Pressure  Trace  for  Calculating  Vdet. 
Occasionally,  speeds  of  30%  to  40%  greater  than  theoretical  velocities  were  seen  in  many 
of  the  experiments.  This  phenomenon  is  explained  by  the  Von-Nuemann  spike  theory. 

He  states  that  there  is  a  short  period  of  over  pressurization  and  overdriven  velocities 
generated  at  the  onset  of  a  detonation  transition. 
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Once  the  wave  speed  is  known  using  this  technique,  it  can  be  compared  to  the 
theoretical  values  to  determine  if  a  detonation  wave  was  present.  The  theoretical 
detonation  wave  velocity  can  be  analytically  determined  using  the  formulas  derived  in 
Chapter  II.C,  however  the  computer  code  TEP  quickly  solves  this  problem  for  a  given 
reactant  mixture,  equivalence  ratio  and  combustor  pressure  and  temperature  parameters 
[Ref.  8].  Figure  [2-13]  is  a  plot  of  this  comparison  completed  for  ethylene/oxygen 
mixture  at  several  equivalence  ratios. 


Figure  [2-13].  Theoretical  vs.  Experimental  Detonation  Velocity  Comparison. 


The  drastic  fall  off  in  wave  velocities  for  the  experimental  cases  at  equivalence  ratios  of 
0.6  and  2.2,  show  the  observed  detonability  limits  of  the  reactants  for  the  geometry 
tested.  The  TEP  code  does  not  account  for  this  limitation  and  carries  out  its  calculations 
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with  no  physical  meaning.  The  error  bars  in  the  plot  quantify  experimental  error  in  the 
time  and  distance  measurement  of  the  propagating  wave  and  vary  with  the  wave's 
velocity.  This  error  analysis  is  rigorously  derived  in  Chapter  IV.B.5.  Once  the 
combustion  event  is  determined  to  be  a  detonation,  the  point  at  which  the  wave 
transitioned  can  be  determined. 

3.  Ignition  Source  Level  and  Location 

The  ignition  source  parameters  are  the  means  by  which  the  combustion  event  is 
initiated  and  play  integral  roles  in  the  development  of  detonations.  The  initial  energy 
provided  is  directly  related  to  how  quick  the  flame  front  becomes  turbulent  resulting  in 
the  deflagration-to-detonation  transition.  This  energy  may  be  provided  by  any  number  of 
means  including  a  spark  ignitor,  focused  acoustic  energy,  or  by  another  detonation  wave. 
The  spark  ignitor  propagates  a  flame  kernel  hemi-spherically  into  the  reactants  initiating 
the  flame  front.  A  detonation  that  sparks  another  detonation  in  a  separate  mixture  is 
called  a  pre-detonator,  and  is  currently  the  means  by  which  the  larger  combustors  of 
PDEs  are  being  initiated.  Each  type  of  ignition  source  drives  a  detonation  transition  at 
different  distances.  The  key  is  to  force  the  detonation  to  transition  as  early  and  simply  as 
possible. 

A  means  by  which  this  is  optimized  is  to  focus  the  available  energy  into  a  smaller 
volume  of  reactants.  This  would  theoretically  create  greater  turbulence  in  the  initial 
flame  front,  driving  the  detonation  to  transition  earlier.  This  research  will  look  into  a 
simple  focusing  technique  of  superimposing  the  flame  fronts  of  an  ignitor  spark.  This 
will  be  accomplished  by  placing  the  ignitor  along  the  combustor's  sidewall  at  different 
locations.  Once  sparked,  the  flame  front  will  propagate  both  forward  and  aft  along  the 
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combustor  until  the  latter  reflects  back  and  coalesces  with  the  forward  wave.  This  creates 
a  greater  energy  flux  through  the  reactants  forcing  a  detonation  to  occur  earlier.  If  proven 
successful,  other  techniques  of  energy  focusing  may  also  prove  to  substantially  decrease 
the  deflagration-to-detonation  transition  distance. 

4.  Thrust  Generated  by  a  Detonation  Wave 

In  order  for  a  detonation  combustion  process  to  be  used  viably  in  a  propulsion 
cycle,  it  must  be  able  to  provide  a  means  of  thrust.  A  look  back  to  the  C-J  cycle  analysis 
and  experimental  results  give  insight  into  this  issue.  A  detonation  in  an  open-ended  tube 
type  combustor  produces  significant  pressure  on  the  combustor's  head  wall.  This 
pressure  is  derived  from  the  post  combustion  properties  of  the  products  behind  the 
detonation  wave  by  performing  a  Taylor-Zeldovich  expansion.  The  properties  of  the  post 
detonation  state  in  the  cycle  can  be  seen  in  Figure  [2-7]  as  state  (2)  and  the  properties  at 
the  head  wall  as  state  (3).  The  pressure  at  the  post  detonation  state  (P2)  in  the  ethylene  / 
oxygen  mixture  is  calculated  to  be  18.7  atm  or  274.9  psia.  The  expanded  pressure  at  the 
head  wall  (P3)  is  calculated  using  the  Taylor-Zeldovich  isentropic  expansion  relation 
shown  by  equation  [2-56].  Using  the  post  combustion  specific  heat  ratio  (72),  P3  is 
calculated  to  be  135  psia. 
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[2-56] 


The  head  wall  pressure  can  also  bee  seen  experimentally  as  shown  by  a  pressure  trace 
three  inches  from  the  head  wall  of  a  detonation  as  it  propagates  through  the  combustor. 
Figure  [2-14]  depicts  this  average  head  wall  pressure  generated  by  the  detonation. 
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Figure  [2-14].  Head  Wall  Pressure  Generated  by  Detonation. 

The  average  pressure  of  the  combustion  products  (P3)  pressure  in  Figure  [2-14],  is  150 
psig.  This  correlates  to  a  momentary  thrust  level  of  1060  lb  for  the  three-inch  diameter 


combustor.  The  gradual  decay  in  pressure  1.35  milliseconds  after  the  detonation  wave 


passes  corresponds  to  the  corresponding  isentropic  blowdown  process  after  the  wave 
exits  the  combustor  into  the  atmosphere.  The  short  time  period  of  the  thrust  is  overcome 
by  the  use  of  a  pulsed  cycle.  A  pulse  detonation  engine  will  need  to  generate  detonations 
at  cycle  rates  of  100  Hz  or  more  to  produce  the  quasi  steady-state  thrust  required  by  a 
flight  vehicle. 
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III.  EXPERIMENTAL  SETUP  AND  PROCEDURE 


A.  INTRODUCTION 

The  test  facility  is  located  at  the  Naval  Postgraduate  School's  Rocket,  Propulsion 
and  Combustion  Laboratory  (RPCL)  in  Monterey,  California.  The  test  apparatus  was 
designed  and  built  specifically  to  investigate  parameters  that  affect  a  transitioning 
detonation  wave.  The  test  cell  in  which  the  experiment  was  set  up  is  a  hardened  bunker 
that  is  well  protected  against  potentially  serious  accidents.  The  facility's  hardware 
consists  of  a  simple  constant  area  tube  combustor,  a  variable  energy  ignition  system, 
control  valves  and  pumps,  and  high-speed  pressure  transducers.  The  gases  that  made  up 
the  reactants  were  a  detonable  mixture  of  ethylene  (C2H4),  oxygen  (O2),  and  nitrogen 
(N2).  The  control  of  all  valves  and  data  acquisition  was  provided  by  a  Microsoft®  Visual 
Basic  GUI  code  developed  specifically  for  this  experiment  [Ref.  9].  The  four 
independent  variables  of  the  experiment,  ignition  energy,  ignition  location,  stoichiometry, 
and  oxygen  (O2)  content,  were  evaluated  in  order  to  determine  their  effect  on  the 
deflagration-to-detonation  transition  distance 

B.  HARDWARE  DESCRIPTION  AND  FUNCTIONS 

1.  Combustor  Tube 

The  most  important  piece  of  hardware  in  the  test  facility  was  the  combustor  tube. 
The  tube  was  75  inches  long  and  3  inches  in  inside  diameter  and  made  of  type-3 16 
stainless  steel.  It  consists  of  two  individual  sections  connected  by  a  flange  assembly  to 
allow  for  different  test  conditions.  Thirty-three  test  ports  are  used  to  place  pressure 
transducers  along  the  inner  surface  of  the  tube  wall.  There  are  also  ports  placed  along  the 
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tube  wall  to  allow  for  injection  and  extraction  of  the  reactants  during  a  test.  The  tube  is 
secured  and  supported  laterally  by  two  notched  braces  and  longitudinally  by  a  thrust  wall 
at  its  head  end  to  absorb  the  impulse  force  generated  by  the  detonation.  Each  brace  is 
mounted  to  a  massive  stationary  test  table.  Figure  [3-1]  is  a  photograph  of  the  combustor 


tube  with  all  its  support  equipment. 


Figure  [3-1 J.  Detonation  Tube  Combustor  with  Support  Equipment. 


2.  Variable  Ignition  System 

An  important  component  evaluated  in  this  experiment  was  the  ability  to 
accurately  vary  the  ignition  source  energy  delivered  to  the  reactants  within  the 
combustor.  This  was  effectively  accomplished  using  the  Unison  manufactured  Vision-8 
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Variable  Ignition  System  (VIS-8).  Figure  [3-2]  is  a  photograph  of  the  Vision-8  system 
with  the  ignitor  installed  into  the  combustor. 


Figure  [3-2].  Unison  Vision-8  Variable  Ignition  System. 


The  ignition  energy  delivered  to  the  ignitor  could  be  varied  from  0.33  to  8.31 
Joules.  The  energy  delivered  to  the  reactants  is  substantially  less  due  to  the  manner  in 
which  a  capacitive  discharge  system  works.  These  losses  are  nominally  75  percent  at 
high  energies  to  65  percent  at  energies  less  than  1  Joule.  The  conditions  for  the  VIS-8 
were  entered  using  a  remote  box  that  allows  the  adjustment  of  energy  levels  and  voltage 
to  achieve  the  desired  ignition  energy. 
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3. 


Vacuum  System 


The  vacuum  system  is  an  integral  part  in  ensuring  the  mixture  being  detonated  is  a 
known  composition.  It  is  made  of  a  solenoid  operated  vacuum  valve  at  the  aft  end  of  the 
combustor  and  the  vacuum  pump  itself.  A  view  of  the  vacuum  valve  looking  towards  the 
head  end  of  the  combustor  can  be  seen  in  Figure  [3-3]. 


Figure  [3-3].  Vacuum  Valve  in  the  Closed  Position. 

The  system  allows  the  combustor  to  be  evacuated  down  to  0.2  psig  so  that  it  may 
be  refilled  with  the  desired  reactants  and  equivalence  ratio.  It  also  allows  the  combustor 
to  be  isolated  from  the  atmosphere  while  the  reactants  are  mixing  through  diffusion  over 
a  10  to  15  minute  period.  Once  the  reactants  were  mixed,  the  vacuum  valve  was  opened. 
A  micro-switch  embedded  in  the  valve  closes  when  the  valve  is  open,  thereby  enabling 
the  ignition  system  to  be  discharged. 
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4.  Fuel  Control  System 

The  fuel  control  system  allowed  the  gases  to  be  delivered  to  the  combustion  tube. 
Since  equivalence  ratio  of  the  reactants  is  a  controlled  variable,  the  measurement  and 
management  of  gas  flow  is  a  critical  factor  in  the  validity  of  the  experimental  procedure. 
It  was  accomplished  through  a  series  of  solenoids,  ball  valves,  flow  restrictive  needle 
valves  and  check  valves.  The  experiment  calls  for  two  separate  pure  gases,  ethylene  and 
oxygen,  and  air,  therefore  two  plumbing  lines  were  run  in  parallel  to  the  combustor. 
Figure  [3-4]  is  a  schematic  of  the  fuel  control  system. 


Figure  [3-4].  Schematic  of  Fuel  Control  System. 


The  gases  were  supplied  from  pressurized  gas  bottles  regulated  down  to  a  constant 
pressure  of  100  psig.  The  combustor  was  isolated  from  the  gases  by  nitrogen-actuated 
ball  valves  that  were  opened  only  during  run  operations.  During  the  filling  process,  the 
specific  gas  was  delivered  at  a  rate  which  resulted  in  a  chamber  pressure  fill  rate  of  0.4 
ps'8/ sec  through  a  needle  valve  into  the  sidewall  of  the  combustor.  Stoichiometry  was 
controlled  by  monitoring  the  tube  pressure.  A  check  valve  provided  the  safety  to  prevent 
backflow  of  gases  or  shock  waves  into  the  fuel  lines.  The  pressures  within  the  lines  and 
in  the  combustor  were  constantly  monitored  through  calibrated  transducers. 
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c. 


REACTANTS 


The  reactants  chosen  for  this  experiments  were  chosen  due  to  the  ease  at  which 
they  detonate  and  the  relatively  short  deflagration  to  detonation  transition  (Xddt) 
distance.  The  mixtures  detonated  were  various  equivalence  ratios  of  ethylene  (C2H4)  and 
oxygen  (O2)  /  nitrogen  (N2)  blends.  The  oxygen  content  of  the  mixture  was  then  lowered 
by  combining  it  with  air  to  test  the  effects  of  nitrogen  dilution  on  Xddt  within  the 
combustor.  The  mixture  was  loaded  into  the  combustor  starting  at  0.3  psi  vacuum  and 
filled  sequentially  until  the  calculated  amount  of  gas  is  present.  Ethylene  was  filled  first 
to  ensure  the  maximum  amount  of  natural  diffusion  could  take  place  under  sub- 
atmospheric  pressures.  This  followed  immediately  by  filling  of  oxygen  until  the  absolute 
pressure  inside  the  combustor  reached  1 .0  atm. 

Using  partial  pressures  of  the  individual  reactants,  derived  from  the  desired 
equivalence  ratio,  the  amount  of  each  reactant  within  the  tube  could  be  accurately 
measured.  The  following  equations  were  derived  for  the  partial  pressure  of  the  reactant 
for  a  given  equivalence  ratio  (cp).  The  equivalence  ratio  is  defined  as  the  ratio  of  the 
actual  fuel  to  oxidizer  ratio  to  the  stoichiometric  fuel  to  oxidizer  ratio.  The  fuel  to 
oxidizer  ratio  (F/O)  can  be  calculated  by  dividing  the  number  of  moles  fuel  by  the 
number  of  moles  of  oxidizer. 
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Using  this  definition  along  with  the  relationship  of  partial  pressures, 
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where  «,■  is  the  number  of  moles  of  the  gas  in  the  mixture,  the  required  partial  pressure  of 
the  reactant  for  a  given  cp  can  be  determined.  For  the  pure  ethylene  /  oxygen  mixture  the 
equation  for  the  fill  pressures  become, 


<p 

Pethyl  <p  +  3Po 


[3-3] 
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The  experiments  investigating  N2  dilution  required  the  use  of  air  within  the  mixture.  The 
significant  percentage  and  known  content  of  nitrogen  (N2)  in  air  lends  itself  as  a  good  gas 
to  incorporate  into  the  reduced  oxygen  mixture.  Experiments  led  to  the  conclusion  that  a 
75/25%  oxygen  to  nitrogen  would  be  the  lowest  possible  ratio  to  evaluate  and  still 
produce  Xddt  distances  within  the  combustor  tube.  Conducting  a  similar  analysis  noted 
above,  the  fill  pressures  for  the  individual  gases  became, 
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Because  the  combustor  is  filled  sequentially,  the  final  pressure  of  the  tube  is  one 
atmosphere  (po).  The  control  screen  that  allows  the  user  to  change  run  conditions  within 
the  program  is  shown  in  Figure  [3-5]. 
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m.  Bun  Conditions 


Set  Ethylene  Partial  Pressure  j:  ■;  ii. 

Air  Partial  Pressure  Jii.C25 
Oxygen  Patial  Pressure  )2.00 

Desired  Equivalence  Ratio  jrCG 

OK  j 

Figure  [3-5].  Control  Screen  for  Altering  Mixture  Ratios.  [From  Ref.  9] 

D.  SOFTWARE  DESCRIPTION  AND  FUNCTIONS 

The  software  used  for  the  experimental  procedure  was  designed  and  written  to 
allow  precise  control  of  the  gas  flow  and  monitoring  of  the  combustor  parameters.  It 
accomplishes  this  while  leaving  the  human  interface  open  to  initiate  the  ignition  or  to 
counter  any  malfunctions  within  the  system  during  the  test.  The  control  code  was 
designed  and  written  using  Microsoft®  Visual  Basic  5.0  GUI  software  [Ref.  9].  Initial 
test  parameters,  including  mixture  type  and  equivalence  ratio  could  be  changed  prior  to 
each  test  within  a  subroutine  of  the  control  code  shown  in  Figure  [3-5].  The  program 
also  allowed  for  the  real-time  visual  representation  of  the  detonation  facility  in  operation. 
All  valve  and  ignition  commands  were  executed  by  TTL  signals  through  a  PI024  board 
with  a  bank  of  Crydom  6321  solid-state  relays  to  the  specific  component.  The  control 
code  continually  updated  the  detonation  tube  pressure  and  the  pressures  of  the  individual 
reactants  to  allow  a  visual  check  of  the  facility's  status.  Imbedded  within  the  program 
were  timing  sequences  that  controlled  critical  safety  measures  in  firing  delays  and  timing 
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of  the  data  acquisition  sequence.  The  layout  of  the  facility  operation  screen  can  be  seen 
in  Figure  [3-6]. 


the  test  setup,  procedure,  data  acquisition,  and  facility  control  prior  to  the  first 
experiment. 

E.  EXPERIMENTAL  TEST  PROCEDURE 

The  testing  procedure  was  a  sequence  of  events  that  allowed  for  careful  control  of 
the  experimental  parameters  while  maintaining  the  physical  checks  required  for  the 
facility's  safe  operation.  In  general,  the  following  is  a  systematic  account  of  the  facility's 
setup,  operation  and  securing  for  a  given  experiment: 
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•  Open  the  facility  and  set  experimental  parameters  including  type  of  mixture, 
equivalence  ratio,  ignition  source  energy  level  and  ignitor  location.  Start  the 
run  sequence  and  allow  code  to  complete  the  following  steps. 

•  Close  the  vacuum  valve  and  turn  on  vacuum  pump  to  decrease  the  tube 
pressure  to  0.3atm. 

•  Open  ethylene  fuel  valves  to  the  combustor  and  fill  to  calculated  partial 
pressure  of  fuel  needed  to  match  the  desired  equivalence  ratio. 

•  Secure  ethylene.  Open  oxygen  lines  to  combustor  and  fill  it  to  a  pressure  of 
latm. 

•  Isolate  the  combustor  and  allow  mixture  to  naturally  diffuse  inside  for  10 
minutes. 

•  Code  enables  ignitor.  Once  a  safety  check  on  the  firing  range  is  cleared,  the 
spark  was  enabled  and  the  vacuum  valve  opened. 

•  The  data  acquisition  board  was  initiated  once  the  end  cap  was  completely 
open. 

•  A  micro-switch  inside  the  vacuum  valve  closes  the  ignitor  loop  when  it  is 
opened,  initiating  the  detonation  with  the  selected  spark  energy. 

•  Upon  completion  of  detonation,  data  was  stored  and  the  facility  closed. 

The  formalized  standard  operating  procedure  for  the  facility's  setup,  operation  and 
securing  for  a  given  test  day  is  included  as  an  Appendix.  Careful  review  of  this 
procedure  prior  to  each  day  of  experimentation  allowed  for  the  safe  operation  and  careful 
control  of  the  testing  integrity. 
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F. 


DATA  GATHERING  AND  PROCESSING 


The  only  source  of  data  that  was  utilized  to  analyze  the  wave  structure  of  the 
combustion  event  was  the  static  pressure  of  the  wave  as  it  passed  over  given  lengths  of 
the  combustor  as  a  function  of  time.  Due  to  the  severe  pressures,  temperatures,  and 
significant  wave  speeds  created  by  a  detonation,  the  pressure  transducers  were  required  to 
be  protected  with  a  silicone  coating.  Seven  Kistler  603B1  pressure  transducers  which 
have  a  frequency  response  greater  than  500  KHz  were  used  to  accomplish  this.  They 
were  amplified  through  501  OB  dual  mode  charge  amplifiers  with  540kHz  notch  filters. 
Two  Microstar  Labs  3400a/415  12-bit  data  acquisition  boards  recorded  the  signals 
generated  by  the  amplifiers.  The  sample  rates  of  the  boards  were  set  to  500  KHz  per 
channel  and  were  synchronized  to  record  all  channels  simultaneously. 

The  data  acquisition  sequence  was  initiated  when  the  opening  vacuum  valve 
closed  the  ignition  loop.  The  board  recorded  data  for  2.5  seconds  before  the  storage 
capacity  on  the  data  acquisition  boards  was  full.  This  allowed  enough  time  to  capture  the 
detonation  wave  as  it  transitioned  and  traveled  out  the  combustor  tube.  Figure  [3-7]  is  a 
plot  of  the  raw  data  points  acquired  through  a  detonation  wave  in  a  typical  detonation 
experiment. 
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Figure  [3-7],  Pressure-Time  Data  Points  Through  a  Detonation  Wave. 


Once  recorded  the  pressure-time  traces  could  be  analyzed  using  the  graphics  software  in 
TecPlot  7.5 

G.  EXPERIMENTAL  TEST  MATRIX 

The  goal  of  the  experiment  was  to  determine  how  specific  changes  in  the  initial 

conditions  of  the  reactants  and  ignition  source  affected  the  length  at  which  the 

combustion  process  transitions  from  a  deflagration  to  a  detonation  wave.  In  order  to 

logically  test  all  parameters  being  investigated  without  convoluting  their  effects,  a  careful 

test  matrix  was  developed.  This  matrix  looked  at  four  independent  variables,  ignition 

source  energy,  ignition  source  location,  reactants’  equivalence  ratio,  and  reactants’ 

oxygen  content.  Table  [3-1]  lists  the  specific  conditions  tested  for  each  parameter. 
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Table  [3-1].  Independent  Parameters  Tested  for  Xppr. 


Variable 

Values  Investigated 

Units 

Ignition  Energy  Level 

0.33, 0.5,  0.8, 1.0, 5.0,  8.31 

Joules 

Equivalence  Ratio 

0.6, 0.8, 1.0, 1.2, 1.5, 2.0 

None 

Ignition  Source  Location 

0.0, 3.0,  7.0 

Inches  from  Head  Wall 

Reactant  O2/N2  Ratio 

100, 75 

Percentage 

Through  combinations  of  each  of  these,  the  deflagration-to-detonation  transition 
(Xddt)  distance  was  determined  by  analyzing  the  developing  wave’s  pressure  traces.  Not 
all  combinations  were  tested,  due  to  time  constraints  and  little  experimental  value  certain 
data  points  would  produce.  Each  experiment  was  conducted  at  least  once,  therefore,  little 
analysis  was  conducted  on  the  shot-to-shot  deviation  of  the  experiments. 


51 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


52 


IV.  EXPERIMENTAL  RESULTS 


A.  ANALYSIS  CRITERIA 

The  object  of  these  experiments  was  to  determine  the  deflagration-to-detonation 
transition  (Xddt)  distance  for  various  combinations  of  initial  conditions  listed  in  Chapter 
III.G.  A  standard  set  of  guidelines  was  outlined  to  compare  experimental  data  for  Xddt 
determination.  The  following  is  a  list  of  these  guidelines  and  an  explanation  of  why  they 
were  positive  indications  of  a  detonation  transition. 

•  The  wave  velocity  is  greater  than  90%  of  the  theoretical  detonation  velocity.  This 
is  the  strongest  indication  that  a  detonation  had  fully  transitioned.  Both  theory 
and  experiments  show  that  there  exists  a  500  to  1200  m/s  difference  in  deflagration 
driven  shock  velocities  and  detonation  wave  velocities.  This  criterion  in  velocity 
difference  is  common  and  has  been  used  in  the  earliest  experiments  of  detonation 
waves  [Ref.  1]. 

•  The  wave's  static  pressure  must  exhibit  a  substantial  step  increase.  The  ZND 
wave  structure  models  the  detonation  wave  with  a  sharp  increase  in  pressure. 

This  is  caused  by  a  normal  shock  followed  immediately  by  a  reaction  zone  as  seen 
in  Figure  [2-3].  Experimentally  this  is  seen  as  a  near  step  increase  in  the  pressure 
traces  of  known  detonations.  This  can  be  seen  in  Figures  [2-8],  [2-9],  [2-1 1]  and 
[3-4],  and  proves  to  be  a  strong  indication  of  a  detonation  wave. 

If  both  of  these  criteria  were  met,  the  wave  was  considered  fully  transitioned  at  that 
transducer  location.  The  Xddt  distance  was  then  deemed  to  occur  midway  between  that 
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transducer  and  the  previous  transducer.  The  margin  of  error  on  the  Xddt  distance 
determination  is  the  distance  between  transducers. 

B.  COMPLETED  EXPERIMENTAL  TEST  MATRICES 
1.  Ethylene  /  Oxygen  with  Ignitor  at  Head  Wall 
The  goal  of  the  first  experimental  test  matrix  was  to  complete  a  detailed  plot  of 
the  effects  of  varying  ignition  energy  level  and  mixture  equivalence  ratio  on  Xddt 
distance.  The  simplest  conditions  were  set  to  ensure  a  detonation  within  the  test 
combustor  and  gain  experience  with  the  gaseous  mixture.  The  reactants  used  were 
ethylene  and  oxygen,  which  are  known  to  have  Xddt  distances  of  less  than  12  inches  near 
stoichiometric  conditions.  The  second  condition  evaluated  was  the  ignition  source 
location.  The  ignitor  was  placed  on  the  head  wall  to  minimize  any  effects  of  multiple 
combustion  fronts  or  shock  focusing.  Pressure  transducers  were  placed  two  inches  apart, 
to  ensure  capturing  estimated  Xddt  distance.  Once  successful  operation  of  the  test 
facility  was  confirmed  and  the  initial  course  test  matrix  completed,  the  test  procedure 
outlined  in  Chapter  III.E  and  the  Appendix  was  initiated. 

The  test  matrix  covered  a  range  of  ignition  energies  from  0.3  to  8.0  Joules  over 
several  mixture  equivalence  ratios  from  0.8  to  2.0.  The  ignition  energy  range  was  based 
on  the  limits  of  the  ignition  system.  The  equivalence  ratio  parameter  range  was  based  on 
the  detonability  limits  found  experimentally  and  discussed  in  Chapter  III.E.2.  The  results 
of  this  experimental  matrix  are  recorded  in  Table  [4-1]. 
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Table  [4-1].  Ethylene  /  Oxygen  Test  Matrix  for  Ignitor  Position  on  Head  Wall 


Energy  Level  Sui 
(Joul 

jplied  to  Ignitor 
es) 

Equivalence  Ratio 

<i> 

Xddt 

1.0 

2.5 

8.0 

O 

oo 

8 

8 

8 

8 

X 

8 

1.0 

6 

6 

5 

5 

5 

5 

1.2 

3 

3 

3 

3 

3 

3 

1.5 

3 

3 

3. 

3 

3 

3 

2.0 

6 

6 

5 

5 

5 

5 

The  data  point  at  a  <|>  of  0.8  and  ignition  energy  of  2.5  was  not  taken,  however  the  point  is 
not  at  an  inflection  and  there  is  no  reason  to  believe  the  trend  would  be  upset.  Each  data 
point  was  only  sampled  once;  therefore,  no  statistical  results  can  be  made. 

Graphically,  this  data  can  be  presented  in  two  separate  ways.  The  first  is  an 
inverted  carpet  plot,  Figure  [4-1],  that  displays  how  Xddt  varies  with  the  initial  mixture 
and  spark  energy  conditions.  The  optimal  conditions  are  also  pointed  out  at  an 
equivalence  ratio  of  1 .3  and  spark  energy  of  0.3  Joules. 
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Figure  [4-1].  Carpet  Plot  of  Ethylene  /  Oxygen  Baseline  Variables  Test. 

The  second  way  to  present  the  data  is  a  contour  plot  of  Xddt  as  a  function  of  the  reactant 
equivalence  ratio  and  ignition  energy.  In  the  plot,  depicted  in  Figure  [4-2],  the  contour 
levels  are  the  transition  distances  for  the  corresponding  ignition  energy  and  equivalence 
ratio.  The  light  color  region  shows  the  optimal  conditions  to  minimize  Xddt  and  required 


ignition  energy. 


Figure  [4-2].  Contour  Plot  of  Ethylene  /  Oxygen  Baseline  Variables  Test. 
Analysis  of  the  results  of  the  experimental  test  matrix  will  be  discussed  in  Sections  C  and 
D  of  this  chapter. 


2.  Ethylene  /  Oxygen-Nitrogen  Blend  with  Ignitor  at  Head  Wall 

The  goal  of  this  test  matrix  was  to  sample  the  effects  of  lowering  the  oxygen 
content  of  the  oxidizer,  through  N2  dilution,  on  XDDt  distance.  The  ignition  level  and 
mixture  equivalence  ratio  range  were  duplicated  and  compared  with  the  first  test  matrix. 
The  ignitor  was  left  on  the  head  wall  of  the  combustor. 
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The  test  matrix  covered  the  same  range  of  ignition  energies  and  mixture 
equivalence  ratios  as  the  ethylene/oxygen  test  matrix.  The  results  of  this  experimental 
matrix  are  recorded  in  Table  [4-2]. 


Table  [4-2].  Ethylene  /  Oxygen-Nitrogen  Blend 
Test  Matrix  for  Ignitor  Position  on  Head  Wall. 


n 

Energy  Level  Suppl 
(Joules] 

ied  to  Ignitor 

Equivalence  Ratio 

4> 

Xddt 

0.33 

0.5 

1.0 

8.31 

0.6 

Outside 

tube 

Outside 

tube 

Outside 

tube 

Outside 

tube 

0.8 

16 

14 

12 

14 

1.0 

16 

12 

12 

10 

1.2 

16 

12 

12 

12 

2.0 

20 

20 

18 

18 

Xddt  transition  distances  could  not  be  determined  for  tests  conducted  at  a  <()  of  0.6. 
Pressure  traces  showed  that  a  strong  shock  wave  was  forming,  but  wave  velocities  at  the 
end  of  the  tube  were  not  yet  equal  to  that  of  a  detonation. 

The  contour  plot  presents  the  results  in  the  best  manner.  The 
are  pointed  out  at  an  equivalence  ratio  of  1.1  and  spark  energy  of  0.5 
is  a  graphical  representation  of  the  contour  plot. 


optimal  conditions 
Joules.  Figure  [4-3] 
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Figure  [4-3].  Contour  Plot  of  Xddt  for  Ethylene  /  Oxygen-Nitrogen  Blend 
with  Ignitor  Located  on  the  Head  Wall. 

3.  Ethylene  /  Oxygen-Nitrogen  Blend  with  Ignitor  3"  from  Head  Wall 

The  goal  of  this  test  matrix  was  to  determine  the  effects  of  relocating  the  ignitor 
on  Xddt  distance.  The  ignitor  was  placed  one  tube  diameter,  3",  from  the  head  wall  to 
determine  the  benifits  of  creating  two  reaction  fronts,  thereby  doubling  the  inital  energy 
release  rate  and  shorten  the  XDdt  distance.  The  ignition  level  and  mixture  equivalence 
ratio  were  again  the  variables  varied  with  the  test  matrix  outlined  in  Section  B.2  of  this 
chapter. 
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The  test  matrix  covers  the  same  range  of  ignition  energies  and  mixture 
equivalence  ratios  as  in  the  previous  experiments.  The  results  of  this  experimental  matrix 
are  recorded  in  Table  [4-3]. 


Table  [4-3].  Ethylene  /  Oxygen-Nitrogen  Blend 
Test  Matrix  for  Ignitor  3"  from  Head  Wall. 


Bi 

Energy  Level  Supplied  to  Ignitor 
(Joules) 

Equivalence  Ratio 

Xddt 

0.33 

0.5 

1.0 

8.31 

0.6 

39 

33.5 

33.5 

33.5 

0.8 

12 

12 

10 

10 

1.0 

12 

10 

8 

8 

1.2 

12 

12 

12 

12 

2.0 

20 

18 

16 

16 

The  contour  plot  of  the  results  is  shown  in  Figure  [4-4].  The  optimal  conditions 
shown  at  the  start  of  the  light  region  are  at  an  equivalence  ratio  of  1 .0  and  spark  energy  of 
1.0  Joules. 
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Reactant  Equivalence  Ratio 

Figure  [4-4].  Contour  Plot  of  Xddt  for  Ethylene  /  Oxygen-Nitrogen  Blend 
with  Ignitor  Located  3"  from  the  Head  Wall. 


4.  Ethylene  /  Oxygen-Nitrogen  Blend  with  Ignitor  7"  from  Head  Wall 


The  goal  of  the  last  test  matrix  was  to  continue  to  evaluate  the  effects  of 


relocating  the  ignitor  on  Xddt  distance.  The  ignitor  was  placed  approximately  two  tube 


diameters,  7",  from  the  head  wall  to  determine  if  there  exists  any  further  benefits  in 


minimizing  Xddt-  The  ignition  level  and  mixture  equivalence  ratio  were  again  varied  as 


described  in  Section  B.2  and  B.3  of  this  chapter. 
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The  test  matrix  covers  the  same  range  of  ignition  energies  and  mixture 
equivalence  ratios  as  in  the  previous  experiments.  The  results  of  this  experimental  setup 
are  recorded  in  Table  [4-4]. 

Table  [4-4].  Ethylene  /  Oxygen-Nitrogen  Blend 


Test  Matrix  for  Ignitor  7"from  Head  Wall. 


H 

Energy  Level  Supplied  to  Ignitor 
(Joules) 

Equivalence  Ratio 

Xddt 

0.33 

0.5 

1.0 

8.31 

0.6 

Outside 

tube 

Outside 

tube 

Outside 

tube 

Outside 

tube 

0.8 

20 

18 

18 

16 

1.0 

18 

16 

16 

16 

1.2 

20 

18 

18 

18 

2.0 

31.5 

31.5 

31.5 

31.5 

The  contour  plot  is  shown  in  Figure  [4-5]  depicts  the  data  contained  in  Table  [4- 
4].  The  optimal  conditions  are  shown  to  be  at  an  equivalence  ratio  of  1.0  and  spark 
energy  of  0.5  Joules. 
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Figure  [4-5].  Contour  Plot  of  Xddt  for  Ethylene  /  Oxygen-Nitrogen  Blend 
with  Ignitor  Located  7"  from  the  Head  Wall. 

5.  Uncertainty  Analysis 

There  are  inherent  errors  in  the  measurements  taken  for  these  experiments  due  to 

the  nature  of  the  parameters  being  measured.  Specifically,  the  two  values  that  are  of 

interest  in  the  uncertainty  calculations  are  in  the  detonation  velocity  (V<jet)  and  the 

deflagration-to-detonation  transition  distance.  Two  parameters  that  are  critical  in  making 

judgments  on  the  detonation  wave  structure  and  transition  characteristics. 

The  detonation  velocity,  Vd*,  was  calculated  by  analyzing  the  pressure-time  trace 

and  applying  the  rate  equation  to  the  passing  detonation  wave.  This  analysis  is  detailed 

63 


in  Chapter  II.E.2.  The  error  in  this  calculation  stems  from  the  two  variables  in  equation 
[4-1].  The  error  in  Ax  is  due  to  the  physical  errors  of  pressure  ports  in  the  combustor. 


[4-1] 


Tolerance  specifications  on  these  ports  were  0.005  inches  from  centerline,  which 
corresponds  to  a  total  distance  error  between  test  ports  of  Ax=±  0.01  inch.  The  second 
source  of  uncertainty  arises  from  the  sample  rate  of  the  data  acquisition  board.  The  board 
samples  at  a  500  KHz  frequency,  therefore  the  pressure  wave  may  hit  the  transducer  in 
between  data  samples.  The  time  error  in  the  rate  equation  becomes  ±0.50x1  O'6  seconds. 
Knowing  these  sources  of  error,  the  detonation  velocity  uncertainty  can  be  calculated 
including  error  terms  from  the  experiments. 

The  error  in  determining  XDdt  is  not  as  definitive  as  the  velocity's  error  term  due 
to  the  required  interpretation  of  the  developing  wave  structure.  This  technique  is 
described  in  detail  in  Chapter  II.E.l.  The  error  in  determining  XDdt  is  derived  from  the 
distance  between  the  pressure  transducers.  The  transitioning  wave  may  not  have  been 
developed  at  one  test  port,  but  at  the  next  show  a  fully  transitioned  wave.  Therefore,  the 
actual  transition  occurred  between  the  two  transducers.  The  error  in  this  determination 
becomes  the  actual  distance  between  the  two  transducers  where  the  transition  occurred. 
This  was  not  a  precise  measurement  with  errors  on  the  order  of  one  or  two  inches  due  to 
human  interpretation  of  the  pressure  traces.  However  trends  can  be  developed  with  these 
measurements  describing  the  effects  of  changing  initial  parameters  on  XDdt 
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C.  EFFECTS  OF  EQUIVALENCE  RATIO  ON  XDDT 


The  detonable  mixture's  equivalence  ratio  was  an  important  factor  in  how  rapidly 
it  was  able  to  transition  into  a  detonation.  This  was  expected  since  most  of  the  energy 
released  during  the  event  comes  from  the  chemistry  of  combustion.  To  analyze  the 
specific  effect  of  altering  the  equivalence  ratio  of  the  reactants,  a  breakdown  of  the 
contour  plots  into  a  two-dimensional  X-Y  plot  is  provided.  By  holding  all  other  variables 
constant,  a  good  indication  of  the  strong  influence  the  mixture  equivalence  ratio  has  on  a 
transitioning  detonation  can  be  seen.  The  following  figures,  Figure  [4-6],  [4-7],  and  [4- 
8],  represent  each  of  the  three  ignitor  locations.  Each  plot  displays  the  four  different 
ignition  energies  tested  for  the  specific  ignition  location.  Five  different  equivalence 
ratios  of  the  ethylene  /  oxygen-nitrogen  blend  were  tested  and  presented. 


F*-  -  E  =  8.31J  E=1 ,0J  —  E=0.5J  •  -  E=0~33J  | 

Figure  [4-6].  Xddt  Distance  as  Function  of  Equivalence  Ratio 
with  Head  Wall  Ignition. 
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|—  •  E=8.31J  •  ■  A  •  •  E=1.0J  —  <—  E-0.5J  •  E=0.33J  | 

Figure  [4-7].  Xddt  Distance  as  Function  of  Equivalence  Ratio 
with  Ignition  3"  from  Head  Wall. 


0.7  0.9  1.1  1.3  1.5  1.7  1.9 

Equivalence  Ratio 

— E=0.33J  ~  m-  E=0.5J  -  -  A -~-~E~=1.0J~—  -  E=8.31J  ] 

Figure  [4-8].  Xddt  Distance  as  Function  of  Equivalence  Ratio 
with  Ignition  7"  from  Head  Wall. 


It  can  be  seen  in  these  plots  that  the  optimized  equivalence  ratio  for  every 
condition  lies  from  a  range  of  cp's  from  0.9  to  1.2.  In  addition,  there  is  a  sharp  increase  in 
Xddt  for  cp's  less  than  0.75.  In  contrast,  there  was  a  much  more  gradual  increase  in  Xddt 
for  higher  equivalence  ratios.  In  a  practical  system,  the  data  supports  the  running  a  Pulse 
Detonation  Engine  with  this  specific  mixture  at  equivalence  ratios  slightly  greater  than 
stoichiometric.  This  would  minimize  the  detonation  transition  distance  while  also 
allowing  for  fuel  load  variations  without  substantially  impacting  the  engine's  cycle. 

These  results  produce  a  repeatable  trend  in  the  behavior  of  a  transitioning  detonation 
wave  varying  mixture  equivalence  ratios. 

D.  EFFECTS  OF  IGNITION  ENERGY  ON  XDDt 

Ignition  energy  was  investigated  to  determine  the  minimum  energy  required  to 
initiate  a  detonation  in  the  shortest  distance.  The  Xddt  distance  was  expected  to  be 
minimized  at  high  ignition  energies  due  to  the  amount  of  energy  initially  supplied  to  the 
reactants.  A  series  of  two-dimensional  figures  show  the  effects  of  ignition  energy  on  a 
transitioning  detonation  for  the  three  geometries  investigated.  The  following  figures, 
Figure  [4-9],  [4-10],  and  [4-1 1],  represent  each  of  the  three  ignitor  locations.  Each  plot 
displays  the  five  different  equivalence  ratios  tested  for  the  specific  ignition  location. 

Four  different  ignition  energies  supplied  to  the  ethylene  /  oxygen-nitrogen  blend  were 
tested  and  plotted.  Ignition  energies  of  2.0  and  5.0  Joules  were  also  tested  and  showed  no 
difference  in  XDDT.  For  this  reason,  data  points  between  1.0  and  8.3 1  Joules  were  not 
shown  as  they  had  little  influence  on  the  final  results. 
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Figure  [4-9].  Xddt  Distance  as  Function  of  Ignition  Energy 
with  Head  Wall  Ignition. 


Figure  [4-10].  Xddt  Distance  as  Function  of  Ignition  Energy 
with  Ignition  3"  from  Head  Wall. 
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Ignition  Energy  Supplied  to  Ignitor  (Joules) 

| phi=0  6  —  -B—  phi=0.8  -  A  -  phi=L0  *  phi=1.2  —  -X  —  phi=2.0  | 

Figure  [4-11].  XDdt  Distance  as  Function  of  Ignition  Energy 
with  Ignition  7"  from  Head  Wall. 

The  data  presented  in  the  above  figures  support  the  fact  that  there  is  little  to  no 

dependence  on  ignition  energies  greater  than  0.5  Joules  on  Xddt-  At  energy  levels  lower 

than  0.5  Joules  the  transition  distance  begins  to  increase  showing  that,  there  is  a 

minimum  ignition  energy  requirement  to  maintain  a  minimum  XDdt  in  the  initiation  of  a 

detonation.  This  result  was  a  common  result  throughout  all  tested  conditions. 

E.  EFFECTS  OF  IGNITION  SOURCE  LOCATION  ON  XDDT 

Three  separate  ignition  locations  were  tested  to  determine  if  there  exists  a  way  to 

increase  the  energy  release  rate  of  the  reactants  in  order  to  shorten  the  deflagration-to 

detonation  transition  distance.  The  first  tests  were  completed  with  the  ignitor  placed  at 

the  head  wall  of  the  combustor.  The  following  two  test  matrices  were  tested  at  ignition 

locations  one  and  two  combustor  tube  diameters  from  the  head  wall.  Reviewing  the  data 

presented  in  Tables  [4-2],  [4-3],  and  [4-4],  it  is  clear  that  there  was  some  influence  of 

ignition  location  on  Xddt-  The  average  value  of  the  minimized  Xddt  for  head  wall 
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ignition  taken  from  Table  [4-2]  is  12.5  inches.  The  average  value  of  the  minimized  Xddt 
for  the  ignition  located  one  tube  diameter  from  the  head  wall,  taken  from  Table  [4-3],  is 
9.5  inches.  The  average  value  of  the  minimized  Xddt  for  the  ignition  located  two  tube 
diameters  from  the  head  wall,  taken  from  Table  [4-4],  is  16.5  inches. 

It  is  clear  that  by  placing  the  ignitor  one  tube  diameter  from  the  head  wall  reduced 
the  Xddt  by  24%.  However,  there  is  a  limit  to  how  far  you  can  move  the  ignitor  from  the 
head  wall.  By  placing  the  ignitor  two  tube  diameters  from  the  head  wall,  the  XDdt  was 
increased  by  24%.  This  result  shows  the  advantage  to  locating  the  ignition  source  to  a 
position  that  amplifies  the  energy  release  rate  of  the  reactants  just  after  ignition.  Future 
testing  could  explore  this  phenomenon  by  altering  ignition  strategies  to  further  reduce  the 
deflagration-to  detonation  transition  distance. 

F.  EFFECTS  OF  DECREASED  OXYGEN  CONTENT  ON  XDDT 

In  order  to  understand  the  effects  of  N2  dilution  in  the  oxidizer,  four  separate 
mixtures  of  ethylene  (C2H4),  oxygen  (O2),  and  nitrogen  (N2)  were  tested  to  determine  the 
effect  on  Xddt-  Since  an  air-breathing  propulsion  system  should  utilize  as  much  air  for 
the  oxidizer  as  possible,  the  pure  oxygen  carried  on  board  must  be  minimized  for  volume, 
safety,  and  performance  considerations.  The  four  mixtures  tested  were  various 
percentages  of  oxygen  and  nitrogen  that  make  up  75%  of  the  total  mixture.  Ethylene  was 
the  fuel  used  for  all  the  tests.  The  first  test  mixture  was  made  up  of  pure  oxygen  with  no 
nitrogen  present.  The  second  test  attempted  to  use  air  (24%  oxygen  to  76%  nitrogen).  The 
third  and  fourth  test  used  a  50%/50%  and  25%/75%  oxygen/nitrogen  blends  respectively. 
Using  these  parameters,  the  effect  of  varying  the  nitrogen  content  in  the  mixture  could  be 
investigated. 
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The  gaseous  mixture  with  no  nitrogen  dilution,  only  ethylene  and  oxygen, 
produced  the  shortest  Xddt,  as  expected.  Averaging  the  minimized  distances  from  Table 
[4-1],  the  tests  show  an  average  XDdt  of  3  inches.  Tests  using  ethylene  and  air  and  the 
50%/50%  oxygen/nitrogen  blend  were  unsuccessful  in  producing  a  detonation  within  the 
combustor  tube.  It  can  be  deduced  from  these  failures  that  if  the  mixture  was  detonable, 
the  transition  distance  is  greater  than  75  inches  and  not  likely  practical  for  the  pre¬ 
detonator  component  of  a  propulsion  system.  The  fourth  and  final  oxidizer  ratio  of 
75%/25%  oxygen/nitrogen  was  tested  and  detonated  within  the  combustor  tube. 
Averaging  the  minimized  distances  from  Table  [4-2],  the  tests  show  an  average  XDdt  of 
12.5  inches.  These  results  show  a  substantial  and  detrimental  effect  of  adding  an  inert 
gas  to  the  mixture  on  its  detonation  characteristics.  Although  not  unexpected,  they  do 
show  the  need  for  alternate  tube  geometries  and  ignition  strategies  to  further  minimize 
Xddt  in  an  air-breathing  Pulse  Detonation  Engine. 
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V.  CONCLUSIONS 


Pulse  detonation  technology  is  theoretically  a  viable  propulsion  concept  with  a 
higher  thermal  efficiency  and  simpler  combustor  design.  Utilizing  detonations  for  an  air- 
breathing  propulsion  system  requires  a  comprehensive  research  effort  that  is  currently 
being  conducted  throughout  the  aerospace  propulsion  industry  and  by  the  Office  of  Naval 
Research.  Pulse  detonation  technology  must  mature  beginning  with  research  in  the 
fundamental  behavior  of  practically  developing  a  detonation  wave  reliably,  many  times  a 
second.  This  thesis  investigated  fundamental  issues  by  exploring  the  effects  of  ignition 
location,  ignition  energy,  and  reactant  mixture  ratios  to  minimize  the  deflagration-to- 
detonation  transition  distance  of  a  single  detonation  wave.  The  test  program  revealed, 
clear  trends  and  demonstrated  the  effects  of  these  parameters. 

The  equivalence  ratio  of  the  ethylene/oxygen/nitrogen  blends  was  the  most 
important  variable  since  the  mixture  makeup  delivers  nearly  all  of  the  energy  release  to 
sustain  the  detonation.  It  was  found  that  the  optimized  equivalence  ratio  for  nearly  all 
condition  range  from  cp's  ~  0.9  to  1.2.  There  existed  a  sharp  increase  in  XDDt  for  (p's  less 
than  0.75,  while  proving  less  significant  increases  for  higher  equivalence  ratios.  In  a 
practical  system,  the  data  supports  the  running  a  Pulse  Detonation  Engine  with  mixtures 
at  equivalence  ratios  that  are  slightly  fuel  rich. 

The  effects  of  ignition  energy  and  ignition  location  on  XDDt  were  conducted  to 
explore  the  minimum  ignition  requirements  to  initiate  a  detonation  in  a  reasonable 
distance.  The  ignition  energy  supplied  to  the  ignitor  was  independently  varied  from  0.33 
to  8.31  Joules  and  revealed  that  ignition  energies  above  0.5  Joules  had  little  to  no  affect 
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on  Xddt-  Only  for  values  below  0.5  Joules  did  the  transition  distance  begin  to  increase. 
This  indicates  that  large  power  requirements  for  the  ignition  system  in  a  Pulse  Detonation 
Engine  are  not  required. 

The  second  part  of  this  thesis  investigated  ignition  locations  and  revealed  that 
when  the  ignitor  is  strategically  placed  one  tube  diameter  from  the  head  wall,  the 
transition  distance  could  be  furthered  reduced.  These  results  lead  to  the  conclusion  that  a 
Pulse  Detonation  Engine  operating  with  the  tested  reactants  could  be  initiated  from 
energy  levels  on  the  order  of  standard  automobile  spark  plugs. 

The  third  and  final  goal  of  this  thesis  was  to  determine  the  effect  of  lowering  the 
oxygen  content  of  the  reactant  on  Xddt.  The  purpose  of  this  investigation  was  to 
determine  the  oxygen  content  required  to  develop  a  detonation  in  a  reasonable  combustor 
length  through  the  deflagration-to-detonation  transition  process.  Four  mixtures  were 
used  with  various  percentages  of  oxygen  and  nitrogen  that  make  up  the  oxidizer  in  the 
reactants.  It  was  found  that  for  the  given  combustor  length  of  75  inches,  the  minimum 
oxygen  to  nitrogen  ratio  was  75%.  This  result  reveals  the  sensitivity  of  detonation 
initiation  in  fuel/air  mixtures  and  must  be  further  studied  to  understand  techniques  which 
may  be  used  to  accelerate  the  process  on  an  actual  PDE. 

This  research  bounds  some  of  the  fundamental  limitations  on  the  initiation  of  a 
detonation  in  ethylene/oxygen/nitrogen  gaseous  mixtures.  Though  not  immediately 
practical  for  real-world  applications,  this  fundamental  work  will  aid  in  future  designs  of 
pre-detonators  for  larger  complex  Pulse  Detonation  Engines. 
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APPENDIX,  DETONATOR  TEST  FACILITY  OPERATING  PROCEDURE 


Standard  Operating  Procedure  (SOP) 

Facility  Open  Procedure 

1)  Verify  “Emergency  Shutdown”  switch  is  pushed  IN. 

2)  Turn  on  1 10VAC  power  supply  switch  in  black  wall  mounted  cabinet. 

3)  Turn  on  yellow  “compressed  gas  present”  warning  lights. 

4)  OPEN  C2H4  and  02  isolation  ball  valves. 

5)  Verify  vacuum  line  isolation  ball  valve  is  OPEN. 

6)  Turn  on  Kistler  electronics  and  set  to  proper  gain  and  set  to  OPERATE  mode. 

7)  Turn  on  video  cameras,  and  other  data  recording  as  required  by  experiment. 

8)  Open  all  gas  bottles  (actuator  N2,  Ethylene  QH4  and  02)  and  set  gas  supply  pressures  in  bottle 
room  to:  PN2  >75  psig,  PC2H4  =  90  psig  and  P02  =  90  psig. 

9)  Turn  on  “Boomer”  computer  and  enter  into  DETONATION  Visual  Basic  control  program. 

10)  In  the  program,  enter  the  FACILITY  OPERATION  form. 


Run  Mode 
11) 

12) 

13) 

14) 

15) 

16) 

17) 

18) 

19) 

20) 

21) 


22) 

23) 

24) 

25) 

26) 

27) 

28) 

29) 

30) 

31) 


Verify  “Emergency  Shutdown”  switch  to  the  OUT  position. 

OPEN  the  facility. 

Wait  for  3  minutes  or  more  for  tube  to  vent  any  combustible  mixture. 

Reboot  ignition  system  and  verify  it  displays  default  settings. 

Verify  all  igniter,  solenoid  control,  and  pressure  transducer  lines  are  secure. 

Set  desired  ignition  levels  on  Remote  Ignitor  Box. 

Set  RUN  CONDITIONS  and  ignitor  energy  level  within  the  program. 

Change  FILE  NAME  to  appropriate  serial  number  to  save  data. 

Cycle  End  Cap  and  verify  ignitor  sparks 

Notify  Lab  personnel  that  “Single  Pulse  Detonation  Run  is  about  to  commence”.  Ensure  all 
personnel  return  to  the  Control  Room. 

Ensure  VHS  tapes  are  VCRs.  Start  Facility  VCR  and  set  to  RECORD. 

*****  *  **  *****  *******  WARNING  ******************** 

The  next  step  will  commence  the  run  sequence . 

To  STOP  the  run,  depress  the  ABORT  pushbutton  in  the  program . 

Depress  START  RUN.  Wait  for  gas  mixing  countdown. 

30  seconds  prior  to  detonation,  check  for  golfers  on  video  monitors. 

When  golfers  are  clear,  turn  on  sirens. 

Start  Line  VCR  and  set  to  RECORD. 

Enable  SPARK. 

Press  FIRE  /  ENDCAP  Control  pushbutton. 

Depress  the  SAVE  DATA  pushbutton 
Turn  off  sirens. 

STOP  VCR  recording. 

Upon  completion  of  data  collection,  CLOSE  the  facility. 


Securing  the  Facility 

1)  Depress  “Emergency  Shutdown”  switch  and  verify  it  is  pushed  IN. 

2)  Exit  DETONATION  Visual  Basic  control  program. 

3)  Turn  off  power  to  ignition  system. 

4)  Turn  off  1 10VAC  panel  switch. 

5)  Secure  N2,  C2KU  and  02  gas  bottles.  Turn  OFF  “Compressed  Gas  Present”  warning  lights. 

6)  Close  C2H4  and  02  isolation  ball  valves. 

7)  Turn  off  Kistler  electronics,  cameras  and  VCR. 

8)  Close  door  to  test  cell  #  1 . 
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